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ABSTRACT  OF  DISSERTATION 


TURBULENCE  SPECTRA  IN  THE  BUOYANCY  SUBRANGE 
OF  THERMALLY  STRATIFIED  SHEAR  FLOWS 

A  generalized  eddy-viscosity  approximation  is  used  to  study 
the  turbulence  spectra  of  thermally  stratified  shear  flows.  For  a 


stationary  process  in  the  wave  number  range  investigated--the  buoyancy 
subrange- -under  the  assumption  of  local  homogeneity  of  the  flow,  two 
governing  spectral  equations  with  six  unknowns are, derived, from  the 


equations  of  motion  and  energy., 


In  order  to  reduce  the  number  of  unknowns  to  two  so  that  the 


spectral  equations  can  be  solved,  a  generalized  eddy-viscosity  is 
used  for  expressing  the  integrated  forms  of  the  inertial  transfers 

of  energy  and  temperature  inhomogeneity,  the  shear  stress  and  vertical 


heat  flux  in  terms  of  velocity  spectrum  cj)(k)  and  temperature  spec- 
trum  $  (k) 

Asymptotic  solutions  are  obtained  in  the  buoyancy  subrange 
where  the  local  production  and  local  dissipation  of  turbuleqfr  $pergy 
is  negligible  as  compared  to  the  inertial  transfer  and  vert ical,  heat; 

flux  terms  when  the  flow  conditions  satisfy  the  criterion 


—  /  or 


g  i£Li  «  a  ri.)2 


T  dZ 


e  © 
T 


In  the  buoyancy  subrange  of  stably  stratified  turbulent  flow, 
the  power  law  for  the  velocity  and  temperature  spectra  is  not  univer- 


sal  but  varies  with  the  flow  conditions  in  the  way  (J>(k)  ~  k  and 

m 


ju  11  7 

^(k)  -  k  where  -  —  >_  n  ^3  and  -1  >_  m  —  *  According  to 


the  measurements  of  velocity  spectra  in  the  atmosphere  (Pinus  and 
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Schcherbakova,  1966;  Myrup,  1968) ,  the  dependence  of  the  power  law  on 
the  flow  conditions  was  confirmed.  The  solutions  of  Bolgiano  (1959) 


and  Lumley-Shur  (1964)  are  only  two  particular  cases  of  the  present 
results  under  certain  flow  conditions. 

^  In  the  case  of  th  stratified  turbulent  flow,  the 


velocity  spectrum  exhibits  aLhumr^Ln  the  buoyancy  subrange  as  a  result 


of  the  energy  input  frgifl  the -tpmp^mtnra_fiel^  to  the  velocity  field. 

On  the  left  side  of  this  humpthe  velocity  spectrum  approaches  a  +1 
slope  and  the  temperature- <-pnrtrnm  rhowc—  n — 3—sJ^ipe.  The  measurements 
of  the  velocity  spectra  in  the  atmosphere  (Ivanov  and  Ordanovich,  19b7) 
confirms  this  tendency. 
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Chapter  I 
INTRODUCTION 

Recent  studies  of  locally  isotropic  turbulence  of  homogeneous 
fluids  in  the  inertial  subrange  have  been  helpful  in  providing  solu¬ 
tions  to  several  engineering  problems  which  are  related  to  air  pollu¬ 
tion,  the  long  distance  propagation  of  ultra-high  frequency  radio 
waves  b^  scattering  in  the  ionosphere ,  and  the  safe  structural  desipr) 
of  high  speed  aircraft . 

When  smoke  or  radio-active  material  is  dispersed  by  turbulent 
diffusion,  Tchen  (1959)  has  shown  that  the  dispersion  from  a  point 
source  can  be  related  to  a  function  of  some  power  of  time  for  non- 
stratified  fluids.  Especially,  when  the  -5/3  law  holds  in  the  iner- 

- -  O.  u/S 

tial  subrange  of  locally  isotropic  turbulence,  the  dispersion  of 
particles  is  proportional  to  t3  where  t  indicates  time. 

In  the  second  problem,  if  it  is  assumed  that  scattering  wave 

numbers  fall  in  the  inertial  subrange  of  locally  isotropic  turbulence, 

11/3 

the  scattering  cross  section  exhibits  a  A  dependence  where  X 

is  the  wave  length  of  the  radio  wave  (Bolgiano,  1959^, 

As  to  the  last  problem,  the  vibration  of  aircraft  due  to  the 
atmospheric  turbulence  can  cause  fatigue  of  aircraft  material  and  may 
even  cause  the  aircraft  to  crash  if  the  critical  frequency  of  vibration 
with  respect  to  the  aircraft  is  induced.  Of  course,  a  better  under¬ 
standing  of  the  energy  spectrum  of  atmospheric  turbulence  can  give 
criteria  for  safe  design  of  high  speed  aircraft. 

However,  in  addition  to  the  complexity  of  turbulence,  the 
atmosphere  itself  presents  complications,  i.e.,  the  atmosphere  is 
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usually  thermally  stratified  and  in  a  state  of  shear.  The  turbulence 
changes  its  spectrum  of  energy  in  the  magnitude  or  scale  since  the 
conversion  of  potential  energy  into  or  from  the  kinetic  energy  of  the 
flow  can  increase  or  decrease  the  kinetic  energy  of  turbulence  depend¬ 
ing  upon  whether  the  stratification  is  unstable  or  stable.  Moreover, 
the  spectrum  of  energy  will  change  its  shape  and  form  because  the 
stratification  can  cause  some  anisotropic  effects  on  the  turbulence. 
Thus,  the  turbulence  will  have  directional  properties  and  local  iso¬ 
tropy  can  never  exist  in  the  wave  number  range  where  buoyancy  is  an 
influencing  factor. 

Based  on  the  assumption  that  the  energy  spectrum  depends  only  on 
the  total  dissipation  of  density  fluctuation  by  molecular  effects  in 
the  buoyancy  subrange  of  the  equilibrium  range  of  turbulence,  Bolgiano 
(1959)  reached  a  solution  of  the  energy  spectrum  being  proportional 
to  k  where  k  is  the  wave  number. 

However,  according  to  another  hypothesis,  Lumley  (1964)  obtained 
a  different  spectral  form  in  this  buoyancy  subrange,  since  he  postu¬ 
lated  that  the  energy  spectrum  <j>(k)  and  the  buoyancy  flux  spectrum 

<b  _(k)  are  functions  of  the  local  energy  transfer  flux  e(k)  and  the 
w  I 

local  wave  number  k  and  that  the  spectrum  of  the  buoyancy  flux  in  a 
stably  stratified  flow  is  proportional  to  the  mean  temperature  gra¬ 
dient.  In  this  way,  Lumley  obtained  a  -5  power  law  of  the  wave  number 
in  the  buoyancy  subrange  if  Kolmogorov^  hypothesis  can  be  extended 
to  this  subrange,  i.e.,  energy  spectrum  is  determined  by  k  and  e(k) 

alone,. 

It  is  clear  that  from  the  above  statements,  Bolgiano* s  and 
Lumley* s  results  seem  to  be  mutually  exclusive  at  first  sight.  Hence, 
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the  author's  motivation  will  be  not  only  to  determine  the  discrepancy 
in  their  basic  assumptions  but  to  search  for  the  basic  mechanism  of 
turbulence  in  a  thermally  stratified  turbulent  shear  flow. 
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Chapter  II 
LITERATURE  REVIEW 

In  this  chapter,  previous  studies  of  turbulence  spectra  in  the 
inertial  subrange  are  briefly  reviewed.  Two  hypotheses  given  by 
Bolgiano  and  Lumley-Shur  to  study  the  turbulence  spectra  of  a  stably 
stratified  flow  are  described.  Recent  works  of  Monin,  Gisina,  and 
Pao  are  stated,  and  some  measurements  of  turbulence  spectra  are 
reviewed . 

2 . 1  Locally  Isotropic  Turbulence --Kolmogorov  Hypotheses 

From  the  definition,  turbulence  is  characterized  as  an 
irregular  condition  of  fluid  flow  in  which  fluid  properties  such  as 
vorticity  components  are  distributed  randomly  in  space  and  time. 

Beyond  its  irregularity,  turbulence  as  a  result  of  nonlinear  inter¬ 
action  shows  turbulent  energy  transfer  through  motion  of  the  eddies. 

This  idea  of  turbulent  energy  transfer  is  characterized  by  L.  F. 
Richardson's  rhyme:  "Big  whirls  have  little  whirls,  that  feed  on 
their  velocity,  and  little  whirls  have  lesser  whirls,  and  so  on  to 
viscosity."  This  idea  may  be  stated  in  a  clearer  form  --  the  turbu¬ 
lent  flow  contains  eddies  of  various  sizes  characterized  by  the  non¬ 
linear  interactions  between  eddies.  In  other  words,  the  turbulent 
energy  is  transferred  from  large  eddies  to  smaller  eddies  until  it  is 
dissipated  into  heat  because  of  viscosity. 

In  the  case  of  a  flow  of  high  Reynolds  number,  Kolmogorov  (1941) 
postulated  that  small  eddies  of  turbulence  are  statistically  steady, 
locally  isotropic,  and  independent  of  the  structure  of  large  eddies  of 
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turbulence  from  which  the  small  eddies  are  generated.  He  further 
postulated  that  the  statistical  characteristics  of  small  eddies  of 
locally  isotropic  turbulence  can  uniquely  be  described  by  parameters 
v  the  kinematic  viscosity  and  e  the  total  dissipation  of  turbulent 
energy  by  viscosity. 

In  his  second  hypothesis,  Kolmogorov  postulated  that  in  the 
universal  equilibrium  range  where  small  eddies  lie,  there  exists  a 
subrange  in  which  the  viscosity  effects  are  negligible  and  only  the 
parameter  e  determines  the  turbulence  structure.  Thus,  based  on 
dimensional  arguments,  the  velocity  structure  function,  i.e.,  the 
averaged  square  of  the  difference  of  velocities  at  two  points  sepa- 

~y 

rated  by  a  distance  r  ,  is 

_ _  2  2_ 

|u(x  +  r,  t)  -  u(x,  t)|2  =  a£  e3  r3  ,  <<  r  <<  J lg 

?  (2-1) 
and  equivalently,  the  three  dimensional  energy  spectrum  has  the  form 
in  terms  of  wave  number  k 

2_  _5 

4>(k,t)  =  a  e3  k  3  ,  fc  _1  =  k&  «  k  <<  S^"1  (2.2) 

^  1/4 

in  which  a  and  are  universal  constants,  =  (v^/e)  is 

the  Kolmogorov  length  scale  that  characterizes  a  cut-off  length  scale 
below  which  viscosity  affects  the  turbulence  structure  essentially, 
and  Z  is  the  length  scale  of  energy-containing  eddies.  Histori¬ 
cally,  the  above  stated  universal  function  was  reached  independently 
by  Onsager  (1945,  1949)  and  von  Weizsacker  (1948). 
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Based  on  Kolmogorov* s  hypotheses,  the  turbulent  motion  of 
small  scale  in  the  inertial  subrange  can  be  predicted.  Recent  measure¬ 
ments  of  the  turbulent  energy  spectrum  of  flow  with  high  Reynolds  num¬ 
ber  show  that  the  one-dimensional  energy  spectra  are  proportional  to 


-5/3 


1 


in  the  inertial  subrange.  The  experimentally  evaluated 


for  one-dimensional  spectra  lies  in  the  range  0.48  ±  0.055,  (Pond 

et  al .  ,  1966)  where  is  an  orthogonal  component  of  k  in  the 

1 8 

streamwise  direction  and  =  —  a  derived  from  the  assumption  of 
isotropy.  Examples  of  those  measurements  are  listed  chronologically, 
Gurvich  (1960)  measured  in  a  wind  over  land.  Grant  and  his  colleagues 
(1962a,  1962b)  in  a  tidal  channel.  Pond  and  his  co-workers  (1963, 

1966)  in  wind  over  water  waves,  Gibson  (1963)  in  a  round  jet,  and  Payne 
and  Lumley  (1966)  in  an  atmospheric  surface  layer  by  an  airborne  hot¬ 
wire  anemometer. 

In  case  the  tqpp^rfif  upg  field  is  considered,  Qbukhoff  (1949) 
and  C(ffrsin  (1951)  extended  the  Kolmogorov* s  hypotheses  to  the  tempera¬ 
ture  spectrum,  i.e.,  in  the  inertial  convective  subrange  the  tempera¬ 
ture  spectrum  also  follows  the  -5/3  law  and  has  the  form 


t>TT(k»t)  =  aT  N  e 


5 

3 


k  <<  k  <<  i, 

e  d 


-1 


(2.3) 


where  is  a  universal  constant,  N  is  the  total  dissipation  of 

temperature  fluctuation  by  molecular  transport.  Experimentally, 
Gibson  and  Schwarz  f 19631  showed  the  existence  of  the  -5/3  law  in  the 
inertial  convective  subrange  of  the  temperature  and  concentration 


spectra  measured  behind  grids  in  a  water  tunnel.  Tsvany  MQfifn  also 
found  the  -5/3  law  of  temperature  spectra  in  the  atmospheric  surface 
layer. 
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Although  Eqs.  2.2  and  2.3  give  good  prediction  of  spectra  in 


the  inertial  subrange  of  the  equilibrium  range,  a  detailed  study  of 
the  spectra  covering  the  whole  equilibrium  range  must  take  v  into 
consideration  as  k  approaches  to  and  beyond  k^  .  To  reach  this 
point,  an  additional  assumption  on  the  mechanism  of  the  turbulent 
energy  transfer  must  be  proposed. 


Historically,  Obukhoff  f 194 IT  first  gave  the  assumption  that 


the  energy  transfer  across  the  wave  number  k  is  analogous  to  the 
process  of  the  production  of  turbulent  energy  due  to  the  work  of 
4lQynolds  stress  against  the  mean  motion.  From  the  other  point  of 
ew,  Heisenberg  f 19481  considered  that  the  eddies  with  wave  numbers 
rger  than  k  act  as  turbulent  eddy  viscosity  on  the  eddies  with 

ve  numbers  less  than _ k  .  From  the  eddies  of  wave  number  less 

than  k  energy  is  transferred  to  the  smaller  eddies  with  wave  number 
greater  than  k  .  Assuming  the  local  property  of  the  energy  transfer 
function,  i.e.,  the  energy  transfer  is  only  a  function  of  wave  number 
k  and  the  energy  spectrum  at  this  local  wave  number  k  ,  Kovasznay^ 
(1948)  obtained  some  solution  also.  Using  a  different  approach-- 


cascade  process  approximation,  Pao_  (1965)  obtained  some  solution  for 


the  locally  isotropic  turbulence  at  high  wave  numbers. 

For  the  purpose  of  generalizing  the  problem,  Stewart  and 
Townsend  (1951)  gave  the  assumption  of  generalized  eddy-viscosity 
which  is  actually  expressed  in  the  form  of  a  series.  Due  to  the 
difficulties  involved  in  arriving  at  a  closed  form  for  the  energy 
spectrum  when  a  series  form  of  generalized  eddy-viscosity  approxima¬ 
tion  is  used,  Panchev  (1967)  used  only  one  term  of  this  series  and 


obtained  some  results  for  locally  isotropic  turbulent  flow.  In 
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particular,  one  thing  must  be  noted  that  Heisenberg's  and  Kovasznav's 
approximations  can  be  deduced  from  Panchev's  approach,  and  moreover, 

Pjio's  method  can  also  be  reached  if  some  special  nondimensional  param¬ 
eter  is  introduced.  Panchev's  book  119671  must  be  referred  to  for 
the  details  of  these  relationships. 

Stimulated  by  Panchev’s  work,  the  author  tried  to  extend  the^ 
generalized  eddy-viscosity  approximation  to  the  thermally  stratified 

turbulent  shear  flow  which  will  be  investigated  intensively  in  the 
next  chapter. 

Before  we  study  the  eddy-viscosity  approximation,  some 
limitation  of  this  approximation  hn,  Batchelor  ( 1953) 

objected  that  introduction  of  the  eddy  viscosity  implies  that  the 
smaller  eddies  must  be  statistically  independent  of  the  larger  eddies. 

But  as  k  approaches  to  k^  ,  this  statistical  independence  does 
not  exist  (Hinze ,  1959).  Thus,  at  high  wave  numbers  the  eddy-viscosity 
approximation  cannot  be  valid  as,  on  the  other  hand,  indicated  by  the 
fact  that  -7  law  at  high  wave  numbers  implied  by  the  eddy-viscosity 
approximation  is  unrealistic  because  -7  law  will  mean  the  discontinuity 
of  velocity  derivatives .  However,  according  to  Kolmogorov’s  hypothe¬ 
sis,  this  statistical  independence  mav  be  assumed  in_tho  ipgptial 
subrange  k0  <<  k  <<  ^  ,  and  thus  the  validity  of  ttle  eddv-viscositv 

approximation  will  be  assumed^ 

2 . 2  y  ^  Hypothesq^on  Stably  Stratified 

Turbulent  Flow 

As  described  above,  the  Kolmogorov’s  hypotheses  shed  some  light 
on  turbulent  structure  of  flow  without  any  thermal  effects.  In  the 
atmosphere  the  flow  is  not  only  compressible  but  thermally  stratified. 
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In  case  free  convection  occurs,  the  atmospheric  turbulence  is  excited. 
Even  when  the  atmosphere  is  stably  stratified,  i.e.,  the  lapse  rate  of 
temperature  is  less  than  the  adiabatic  one,  there  is  convincing  evi¬ 
dence,  excepting  when  very  strong  stable  stratification  occurs,  that 
there  exists _a  random,  irregular  mot ion--turbulent  motion  in  the  atmos¬ 
phere  (Kellogg,  1956) .  Due  to  the  existence  of  thermal  stratification, 
the  gravitational  force  must  be  introduced  into  the  equations  of 
motion,  and  the  potential  energy  of  the  flow,  as  an  evidence  of  the 
gravitational  force,  will  affect  the  energy  balance  of  the  flow. 

Hence,  it  can  be  expected  that  any  variation  of  the  turbulent  kinetic 
energy  must  be  a  function  of  atmospheric  thermal  stratification. 


2.2.1  Mechanism  of  turbulence  in  thermally  stratified  flow  - 
As  a  result  of  the  introduction  of  gravitational  effects  due  to  the 
thermal  stratification,  the  turbulent  field  becomes  anisotropic  since 
the  vertical  velocity  fluctuation  is  suppressed  if  the  flow  is  stably 
stratified  and  is  excited  if  the  flow  is  unstably  ^IxaJLi £i It  can 
be  expected  that  in  the  absence  of  shear,  the  turbulent  field  as  a 
first  approximation  tends  to  be  axisymmetric  with  respect  to  the  verti¬ 
cal  axis  to  which  the  gravitational  force  is  oppositely  parallel. 
Incidentally,  the  anisotropic  effect  will  appear  in  the  turbulent 
energy  spectrum  since,  in  the  range  of  wave  numbers  where  damping  or 
excitation  of  turbulence  by  buoyancy  force  occurs,  a  part  of  turbulent 
energy  is  abstracted  from  or  into  the  turbulent  velocity  field  and  is 

converted  into  or  from  potential  energy  depending  on  whether  th^flow 

is  stably  or  unstably  stratified. 


XQ,  lybt xi/Jth  ? 
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For  a  stably  stratified  flow,  energy  drained  out  of  the 
turbulent  velocity  field  may  propagate  away  in  the  form  of  internal 


mot 


•T'  ^ 


gravity  wavf  disturbances  of  Long-Hines  type  (Long.  1953.  1955: 
lines,  1960) may  cause  the  production  of  density  or  temperature 
inhomogeneity  which  is  transferred  inertial ly  to  smaller  eddies  and 
finally  smeared  out  bv  the  nnlprnlar  effects.  Thus,  in  the  wave  num¬ 
ber  range  where  the  buoyancy  force  effects  predominate,  the  turbulent 
pj  a(  f  energy  transfer  c^reas^  with  wave  number  the  transfer  of  den¬ 

sity  or  temperature  inhomogeneities  ^ncreases^  Hereafter,  the  total 
turbulent  energy  dissipation  by  viscosity  is  reduced  and  the 

~  1/4 

Kolmogorov  wave  number  (e/v  )  decreases  accordingly.  In  other 

words,  the  turbulent  scale  at  which  the  viscous  cut-off  fl££LLLS  -wii  1  1 
indirect ly fTncreasey through  the  effects  of  stable,  stratification . 

On  the  other  hand,  in  the  case  of  unstably  stratified. 
turbulent  flow,  the  potential  energy  of  the  temperature  or  density 
field  is  converted  to  the  turbulent  velocity  field  and  the  velocity 
spectrum  may  exhibit  a  hump  in  the  buoyancy  subrange  where  the  gravi¬ 
tational  force  affects  essentially.  Of  course,  the  temperature  or 
density  spectrum  in  this  buoyancy  subrange  may  have  a  steeper  slope 

as  a  result  of  energy  export.  It  can  also  be  expected  that  the  total 
dissipation  of  energy  increases  and  the  wave  number  a±. 
of  the  molecular  effects  occurs  is  increased. 


Keeping  the  above  described  mechanism  of  stratified  turbulent 
flow  in  mind,  Bolgiano's  and  Lumley-Shur f s  hypotheses  on  the  stably 
stratified  flow  are  now  introduced. 
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2.2.2  Bolgiano's  hypothesis  -  If  the  Reynolds  number  of  a 
stably  stratified  flow  is  sufficiently  large,,  Bolgiany  1 19591  postu¬ 
lated  that  the  equilibrium  range  of  the  turbulent  energy  spectrum  can 
be  divided  into  three  distinct  subranges: 

(1)  the  buoyancy  subrange  in  which  the  larger,  anisotropic 
eddies  are  directly  influenced  by  the  density  stratifi¬ 
cation  . 

& i  the  inertial  subrange  in  which  the  anisotropic  effects 

due  to  buoyancy  force  decrease  rapidly  and  the  classical 
hypothesis  of  locally  isotropic  turbulence  is  applicable, 
and 

(3)  the  difssipatiop  subrange  at  high  wave  numbers  where  the 
molecular  effects  dominate. 

It  is  obvious  that  the  last  two  subranges  fall  into_the  cate- 
gory  of  the  locally  isotropic  turbulence.  But,  the  buoyancy  subrange 
needs  special  analysis  and  great  attention.  As  a  lower  limit  of 
scale  cut-off  for  the  buoyancy  subrange ,^BoTgiano^introduced  the 
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where  N  is  the  total  dissipation  of  temperature  fluctuation  by 
molecular  transport,  is  the  total  turbulent  energy  dissipation 

by  viscosity,  and  J3  =  g/T  in  which  g  is  the  acceleration  of 
gravity  and  T  the  mean  temperature. 

Thus,  Bolgiano  further  postulated  that  there  exists  a  wide 
range  of  wave  number  between  the  scale  of  energy-containing  eddies 
and  the  Obukhoff  length  ,  or  equivalently,  there  exists  a 
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buoyancy  subrange.  He  also  assumed  that  in  this  wave  number  range  je^ 
is  comparatively  much  smaller  than  e(k)  the  local  rate  of  the  iner¬ 
tial  transfer  of  turbulent  energy  (energy  transfer  flux)  and  that  the 
local  dissipation  in  this  subrange  is  so  small  that  the  statistical 
properties  of  turbulence  such  as  velocity  and  temperature  spectra  are 
only  a  function  of  N  ,  6  and  wave  number  k  .  Dimensional  argument 
gives  readily 

k  «  k  «  l  _1  .  (2.4) 

e  o 

y  transfer  flux  e(k)  approaches  to 
-5/3  law  holds 

<t>(k)  ~  e3  k  3  ,  i l  ~1  «  k  <<  £d-1  .  (2.5) 


Similarly  the  temperature  spectrum  can  be  worked  out  dimensionally  as 
a  form 


k  <<  k  << 
e 


(2.6) 


4  _1  <<  k  «  St'1  .  (2.7) 

o  d 

2.2.3  Lumley-Shur ' s  hypothesis  -  In  contrast  to  Bolgiano's 
hypothesis,  Lumley  (1964,  1965)  developed  a  new  hypothesis  for  the 
turbulence  spectrum  of  a  stably  stratified  flow.  These  two  theories 
are  mutually  exclusive  since  they  are  based  on  entirely  different 
physical  backgrounds  and,  of  course,  lead  to  different  predictions 


and 


<t>TT(k) 


N  e 
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for  the  spectral  forms  of  turbulent  energy  and  temperature  fluctuations 
in  the  "buoyancy"  or  "locally  inertial"  subrange. 

In  his  paper,  Lumley  first  extended  the  original  Kolmogorov 
hypotheses  in  the  inertial  subrange  into  the  "locally  inertial"  sub- 

range.  From  Kolmogorov’s  hypothesis,  the  statistical  properties  of 
turbulence  in  the  inertial  subrange,  if  the  inertial  subrange  exists, 
are  characterized  uniquely  by  e  and  k  as  stated  in  section  2.1. 

In  the  inertial  subrange,  e  ,  the  turbulent  energy  dissipation  by 
viscosity,  is  in  fact  the  energy  transfer  flux  through  wave  numbers. 
Thus ,  Lumley  postulated  that  the  statistical  properties  of  turbulence 
of  a  stably  stratified  flow  in  the  inertial -buoyancy  subrange,  such 
as  energy  spectrum  <Hk)  and  heat  flux  spectrum  ^wy(k)  >  are  deter- 
mined  by  the  wave  number  k_  and  the  local  energy  transfer  flux  e(k) 
at  this  wave  number  k 

In  addition  to  the  above  hypothesis  as  an  extension  of 
Kolmogorov’s  hypothesis,  Lumley  postulated  further  that  the  temper a- 
ture  fluctuation  field  is  determined  solely  by  the  velocity  field. 

From  the  temperature  fluctuation  equation  of  flow  with  high  Reynolds 
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in  which  x^  are  the  Cartesian  coordinates,  x^  is  in  the  streamwise 


direction,  x0  lateral  and  x^  vertical;  0  is  the  fluctuating 
temperature;  LL  are  the  instantaneous  velocities ;  T  is  the  mean 
temperature;  the  repeated  index  denotes  summation  and  L^  is  the 
Eulerian  space  integral  scale,  Lumley  derived  that 
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where 

reach 


(2.9) 


a(x,  t)  is  the  position  at  t=0  of  a  particle  which  will 
x  at  time  t  .  Thus, 


completely. 

Now,  for  simplicity,  assume  the  mean  temperature  gradient 
exists  only  in  the  vertical  direction,  i.e.,  8T/8x^  =  3T/3x ^  =  0 

Then  the  heat  flux  has  only  the  form  0Uj  'ID  4>w^(k)  the  spectrum 
of  the  vertical  heat  flux  is  integrated  spherically,  the  direc¬ 

tional  information  can  be  missed  for  simplifying  derivations.  Thus, 


in  which  ^(k)  is  the  spectrum  form  of^A^j,do  is  the  surface 

element  of  a  sphere  with  radius  k  =  |£|  ,  and  dT/dz  =  dT/dx 


Finally,  based  on  the  above  derivation,  Lumley  postulated 
that  the  spherically  averaged  spectrum  of  the  vertical  heat  flux 
<frwT(k)  is  proportional  to  the  mean  temperature  gradient  in  vertical 
direction  dT/dz  for  a  stably  stratified  flow.  ©  from  the  dimen¬ 
sional  reasoning,, 
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(2.11) 


in  which  a  is  constant. 

Now,  if  the  production  of  turbulence  is  very  weak  in  the  wave 
number  range  considered,  we  have 


4-  not 
f  cR 

K*  •  , 

ft/**?"' 


(2.12) 


which  results  in 

2  2 
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(2.13) 


where 


_g_  £L  P  3 

T  dz 

i 


2  1 
4 


(2.14) 


After  inserting  e(k)  into  the  generalized  Kolmogorov's 


spectrum 


<f>(k)  =  ae3  (k) k  3 


k  <<  k  <<  a. 

e  d 


(2.15) 


-1 


N 


(2.16) 
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Based  on  the  same  argument  described  above,  Luirtlev^Shur^^ 


(2.17) 


7dT  2 
eT  -T~ 
dz 


gN 


«  1 


In  order  to  make  the 


(2.18) 


Now,  the  discrepancies  between  Bolgiano's  and  Lumley-Shur ' s 
hypotheses  are  apparent  when  Eqs.  2.4,  2.6,  2.16  and  2.17  are  reviewed. 
Their  differences  in  physical  background  can  readily  be  seen  from 
their  hypotheses.  The  details  will  be  discussed  and  compared  after 
the  generalized  eddy-viscosity  approximation  is  introduced  in 
Chapter  III. 


2 . 3  Mon ini s .  Gisina's.  and  Pao's  Works 

Besides  the  above  stated  Bolgiano1 s  and  Lumlev-Shur's  hypotheses, 
Mpn in  (1962)  and  (^i^ina  (1966)  used  Heisenberg's  eddv-vi srosi tv. 
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approximation  to  study  the  turbulence  spectrum  in  a  stably  stratified 


flow.  Moninjobtained  Bolgiano's  solution  in  the  buoyancy  subrange 
for  the  stably  stratified  flow,  and  found  some  humps  in  the  spectra 
in  the  case  of  unstable  stratification.  Gisind  considered  more  com¬ 
plicated  conditions  of  stably  stratified  flow: 

(1)  weak  interaction  of  velocity  and  temperature  fields, 

(2)  strong  interaction  of  velocity  fields  and  temperature 
fields , 

(3)  strong  interaction  of  temperature  fields  and  weak  inter¬ 
action  of  velocity  fields . 


The  results  for  Gisina's  first  two  conditions  are  trivial 
since  Tchen's  arguments  f!953)  on  the  shear  flow  can  be  applied  to 
the  thermally  stratified  shear  flow,  and  as  can  be  expected,  the 


/ 


spectra  of  velocity  and  temperature  are  proportional  to 

iitiorJi  and  a  -1  law  i 


•5/3 


|Hk  ~  for 


ther first  flow  conditi 


s  obtained  for  both  snectra. 


,  dT 

>>  b  £  ~T~ 
dz 


of  velocity  and  temperature  in  case  of  the  (secondly  1 pw  condi t  .ion.  As 
to  the  last  flow  condition.  Gisina  obtained  Bolgiano's  solution  in 
case  the  relationship  between  parameters 


(2.19) 


in  which  b  is  a  numerical  constant  corresponding  to  the  ratio  of^ 
kinematic  eddy  viscosity  of  momentum  to  kinematic  .eddy  conductivity 
of  heat  flux^  can  be  fulfilled.  As  we  can  see  later.  Mopin'  s  and 
Gisina's  solutions  are  only  some  special  conditions  of  the  generalized 
eddy-viscosity  approximation  considered  injthis  present  study,  thus , 
their  works  will  not  be  reviewed  in  detail. 
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In  Pao's  na:)dr  ( }9b7)  ~the  cascade  process  is  applied  to  solve 
a  thermally  stratified  shear  flow,  however,  he  did  not  obtain  any 
power  law  in  the  buoyancy  subrange.  Now,  the  situation  is  very  clear, 
different  approaches  used  by  different  people  to  study  a  problem--a 
thermally  stratified  turbulent  shear  flow  problem  result  in  different 
solutions.  If  they  do  not  contradict  one  another,  there  must  exist 
some  way  to  solve  this  problem  and  to  explain  the  discrepancies  among 
them.  In  the  following,  the  author  investigates  this  problem  by  means 
of  the  generalized  eddy- viscosity  approximation.  Before  reaching  this 
point,  some  turbulence  spectral  measurements  in  the  atmosphere  and  in 
a  wind  tunnel  will  be  reviewed. 

2.4  fjfia^urements  of  T^rfcyl^p^e  Spcotr/i  in  the  Atpiospb^je 

Turbulence  spectra  have  been  measured  in  the  surface  layer 
and  in  the  free  atmosphere  by  several  authors.  Different  stratifica¬ 
tions  of  flows  were  involved  in  these  measurements.  In  case  of  neu- 
tral  stratification,  the  -  5/3  law  holds  for  a  wide  range  of  wave 
numbers  as  can  be  expected  from  Kolmogorov  Vs^hypotheses  for  locally 
isotropic  turbulent  flow  of  neutrally  stratified  fluids . 

Figure  la  displays  a  spectral  density  curve  of  longitudinal 
velocity  component  taken  at  500  m  above  the  ground  when  the  lapse 
rate  of  temperature  from  ground  to  1000  m  is  1°C  per  100  m,  i.e.,  the 
adiabatic  lapse  rate,  and  the  mean  velocity  gradient  is  0.36  m  per  sec 
per  100  m  (Pinus  and  Shcherbakova,  1966) . 

In  the  surface  layer  the  -  5/3  law  was  also  observed  for  the 
energy  spectrum,  the  lower  limit  of  the  inertial  subrange  is 

related  to  stratification  or  Richardson  numbe:y .  Generally  speaking. 
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the  stable  stratification  shrinks  the  inertial  subrange  and  the 

unstable  stratification  extends  it  (Zubkovskii,  1962,  Gurvich ,  JL960] .  ^ - 

SET 

In  the  free  atmosphere  the  spectra  has  a  more  complicated  ■ _ ■ 

form  because  of  buovancy  effects .  In  the  following,  two  categories 
of  flows  are  described--stab le  and  unstable  strat if icat ions--but  the 
details  of  the  discussion  will  be  in  Chapter  IV  where  the  numerical 
solutions  of  the  present  study  are  given. 

2.4.1  Stable  stratification  -  As  described  in  section  2.2.1, 
some  energy  will  be  abstracted  from  the  velocity  field  and  fed  into  the 
temperature  or  density  field.  It  can  be  expected  that  the  velocity 
spectrum  in  the  buoyancy  subrange  will  present  a  steeper  slope  than 
-  5/3  .  Shur  (1962)  first  showed  the  existence  of  the  buoyancy  sub¬ 
range  from  his  measurements.  Later,  Tp*nus  and  Shcherbakova  (1966] 
measured  the  velocity  spectrum  in  the  atmospheric  layer  from  400  to 
4500  m.  In  case  of  stable  stratification,  the  slope  of  measured 
velocity  spectra  in  the  buoyancy  subrange  increased  with  height  for 
roughly  the  same  mean  temperature  gradient.  The  exponent  n  of  the 
velocity  spectrum  k  n  in  the  buoyancy  subrange  varies  from  2.0  to 
3.5.  For  the  sake  of  interest,  the  following  Table  1  digested  from 
Pinus  and  Shcherbakova  ( 1966)  is  listed.  For  better  understanding, 
three  typical  energy  spectra  of  stably  stratified  flows  from  their 
measurements  are  demonstrated  in  Figs,  lc.  Id  and  le. 

In  1963,  Pinus  measured  the  spectral  density  of  the  horizontal 
velocity  component  at  heights  of  6-12  km.  Fig.  2b  displays  one  plot 
of  his  results  which  appears  steeper  slope  than  -  5/3  in  a  certain 
wave  number  range, of  suectyum.  Another  measurement  by  Vinnichenko 
(1966)  as  displayed  in  Fig.  2a  indicates  the  existence  of  the  buoyancy 
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subrange.  Recently,  Myrup  (1968)  found  the  buoyancy  subrange  in  his 
measurements  due  to  the  fact  that  the  steepened  slone  is  close  to  -3 
for  the  longitudinal  velocity  fluctuations  and  between  -  11/5  and  -3 
for  the  vertical  velocity  fluctuation. 


TABLE  1 


Height , 

m 

Vertical 
Temperature 
Gradient , 

deg/100  m 

Vertical  Gradient 
of  the  Mean  Wind 
Velocity, 

m  sec“Vl00  m 

Number 

of 

Spectra 

Range  of 
Variation 
of  n 

n 

400-700 

0.65 

1.84 

9 

2. 0-2. 9 

2.43 

700-1200 

0.61 

0.78 

17 

2. 0-2. 9 

2.50 

1200-1700 

0.76 

0.86 

6 

2. 2-3.5 

2.83 

1700-2500 

0.46 

0.45 

6 

2. 3-3. 3 

2.70 

2500-3500 

— 

— 

1 

2.8 

— 

3500-4500 

0.72 

1.13 

9 

2. 7-3. 5 

3.10 

2. 

.4.2  Unstable 

stratification  -  In 

this  case,  the  velocity 

field  absorbs  energy  from  the  temperature  field  as  potential  energy 
is  converted  into  kinetic  energy;  some  hump  in  the  velocity  spectrum 
can  be  expected.  Fig.  lb  shows  a  spectral  density  curve  of  longitudi¬ 
nal  velocity  component.  (Pinus  and  Shcherbakova,  1966).  In  1967, 

Ivanov  and  Ordanovich  made  a  more  detailed  investigation  of  velocity 
spectra  for  unstable  stratification  in  the  low  frequency  range.  In 
Fig.  3,  some  typical  examples  of  the  measured  velocity  spectra  are 
presented.  (Ivanov  and  Ordanovich,  1967). 

As  to  the  temperature  spectrum,  lesser  information  is  available. 
Tsvang  (  3  96^)  measured  some  spectral  density  curves  of  temperature 
for  both  stable  and  unstable  stratifications.  From  his  measurements, 
deviations  from  -  5/3  can  be  seen  as  the  wave  number  is  less  than 
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-4  -I 

about  10  cm  for  either  stratification,  however,  no  detailed 
discussion  on  these  deviations  from  the  theoretical  point  of  view  has 
been  given  in  his  paper. 

2 . 5  Measurements  of  Turbulence  Spectra  in 

Although  the  inertial  subrange  for  locally  isotropic  turbulence 
has  been  confirmed  as  described  in  section  2.1,  the  buoyancy  subrange 
has  not  been  obtained  in  the  laboratory  since  it  is  still  difficult 
to  generate  turbulence  of  laboratory  scale  with  buoyancy  effects. 

Cermak  and  Chuanv  f 19651  measured  some  vertical  velocity  spectra  in 
thermally  stratified  shear  flows,  however,  no  buoyancy  subrange  was 
observed  since  the  ’’buoyancy  subrange”  mentioned  in  their  paper  lies  in 
the  viscous  dissipation  subrange  evidently.  Also,  Arya  f!9681_did  not 
find  any  buoyancy  subrange;  however,  the  vertical  velocity  and  tempera¬ 
ture  spectra,  measured  at  close  wall  regions  where  both  velocity  and 
temperature  gradients  are  great,  present  -1  slope  at  lower  wave  number 
range  as  predicted  by  Gisina  (1966). 
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Chapter  III 
THEORETICAL  STUDY 


Most  of  the  previous  studies  reviewed  were  made  to  find  some 
asymptotic.  solu^ion^  for  which  some  specific  restrictions  are  assigned 
to  the  flow  conditions.  Thus,  it  would  be  helpful  to  retain  every 
factor  in  the  energy  balance  equation  and  to  find  some  continuous 
spectra  if  the  whole  view  of  turbulence  structure  is  to  be  obtained. 
rIn  this  chapter,  the  spectral  equations  of  the  turbulent  energy  and 
the  temperature  fluctuation  are  derived;  the  generalized  eddy-viscosity 
approximation  is  introduced ;  and  some  analytical  and  asymptotic  solu- 
tions  will  be  given. 


3.1  Deyivat^op  of  the  SjuapfrrftJ  Equations  of  Turbulence  Energy  and 
Temperft^ire  Fluctuation 

In  an  incompressible  turbulent  flow,  the  Navier-Stokes 
equation  is  assumed  to  be  the  equation  governing  the  variation  of  the 
spatial  distribution  of  the  velocity  with  time.  For  simplifying  the 
derivation,  dynamic  v 1  s cysity.. $$ sumed  to  be  cons t  an  t  and  Boussinesq  Ts 


where  IL  is  the  velocity  vector, 
mean  density  of  the  flow  field,  P 


P 


is 


is  the  density,-^  pq 
the  pressure,  and  \i 


*  (3.2) 

is  the 
is  the 
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dynamic  viscosity.  The  assumption  for  the  incompressibility  of  the 


t u rbulent  flow  js  unnecessarv^slnce  ^]ie  density  can  be  replaced  by  the 
jotential  density  if  a  compressible  fluid  like  atmosphere  is  concerned 
In  that  case,  we  need  only  assume  the  flow  speed  is  small  compared  to 
the  speed  of  sound,  or  equivalently  to  say  that  the  Mach  number  is 
much  less  than  1.  The  derailed  study  can  be  referred  to  the  papers 
by  Long  (1953b) ,  Batchelor  (1953b),  Bolgiano  (1962),  Lumley  and 
Panofsky  (1964).  However,  for  the  sake  of  simplicity  in  our  deriva¬ 
tion,  incompressibility  is  assumed  for  the  turbulent  flow. 

In  addition  to  the  governing  equations  for  the  velocity  field, 
an  equation  for  the  temperature  field  which  causes  the  density  flue- 


A 


/ 


f  9T  T1  91 

9 

/  3T  1 

1 

(3.3) 

Po  Cp  l  9t  i  9x. 

UT  3x. 

3x. 

f 

\  i  i 

-T  1 

1  i  1 

J 

where  c^  is  the  specific  heat  capacity  at  constant  pressure,  T  is 

the  temperature  and  u  is  the  thermal  conductivity,  y  and  c 

1  1  *  P 

are  assumed  to  be  constant. 

For  a  turbulent  flow,  the  fluid  properties  can  be  split  into 
two  parts: 


p/FFUS/OX 

Ci.t, 

o?  re  hi 


u.  =  u.  +  u. 

iii 


p  =  p  +  p 


i 


,  p  =  p  +  Pj  , 


T  =  T  +  0 


where  the  bar  denotes  the  t ime^average  or  ensemble  average,  and  u^  , 
p  ,  p^  ,  and  0  are  the  fluctuations  about  their  corresponding 
averages . 
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/  9U .  3u .  _  3U.  31) .  _  3u.  3u. 

Po \TT  +  3t^  +  Uj  3x~  +  Uj  3x7  +  uj  3^T  +  uj  33T 
'  3  3  3  3 


3  CP  +  Px) 


3x. 

i 


-  9  l 

-  (p  +  pjg,  +  u  ■  +  y- 


32u. 


r6i  M  3x.3x.  3x .  3x . 

3  3  3  3 


3U . 


3u. 
+  — L 

dx.  3x. 

1  1 


(3.4) 

(3.5) 


and 


P  c 
O  p 


3T  30  rr  3T  3T  —  30  30 

-r—  +  77—  +  U.  - —  +  u.  - —  +  U.  -  +  u 


3t  3t  i  3x.  i  3x.  i  3x.  i  3x. 

l  l  l  l  / 


=  u 


32T 


+  y  n 


320 


T  3x.3x.  T  3x  .  3x . 

J  3  3  3 


(3.6) 


Assuming  the  mean  flow  is  stationary  and  taking  the  average 
(either  ensemble  or  time  average)  of  Eqs.  3.4,  3.5,  and  3.6,  we  have 
thus 

,  a77  ?hi  _ 

3P  - 

Si  '  »-> . 

3  3 


313 .  3u. 

77  i  1 

U  .  - —  +  u  . 


j  3x.  j  3x . 

3  3  ' 


„  a2ui 

dx.  h  ^i  ^  dx .dx. 


(3.7) 


dU. 

i 

dx . 

l 


=  0 


(3.8) 


P  c 
O  p 


77  3T  30 

U.  - —  +  u. 


i  3x. 

l 


i  3x. 

l 


=  U 


32T 


T  3x . 3x . 
J  J 


(3.9) 


Jjuj^y^ac^j^g  Eqs .  from  Eqs.  ^ . 4 ,  3.5  and ^3^6^  respec¬ 

tively,  the  turbulent  equations  for  velocity  and  temperature  fields 


become 
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3u. 


9U. 

3u. 

3u. 

3<*U3 

l 
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j 

J 

9P1 

gi  + 
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i 

ax.  p  l 
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1 

j  j 
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3u. 

1 


j  / 


3u. 
_ 1 

3x. 

i 


,30  3T  77  30  30 
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o  p  3t  i  3x.  j  3x .  j  3x. 
'  J  J  J 


-  u. 


30 


i  SX.  , 


=  u 


a2e 


T  3x . 3x . 
J  J 


(3.10) 


(3.11) 


(3.12) 


Taking  the  advantage  of  incompressibility,  Eqs .  3.10,  3.11, 
and  3.12  are  rewritten  as 


in  which  v  =  y/p  is  the  kinematic  viscosity,  v 


u„/p  c  is 
T  o  p 


the  thermal  diffus 


ivity ,  and  0/Tq 

Vmmb 


-  p./p  comes  from  the  as 

1  o  i  - 


sump¬ 


tions  that  the  flow  is  incompressible  and  the  (c^-cv)pT 


holds,  where  cy  is  specific  heat  capacity  at  constant  volume. 
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Multiplying  Eqs.  3.13  and  3.14  by  u'^  and  0^  respectively, 
where  the  prime  denotes  that  the  quantities  are  measured  at  x'=x+r  , 


(3.15) 

BMBACi 

■f 

pUfCT- 


(3.16) 


multiplied  by  u^  and  0  at  x  ,  then  we  have 


3u' 


9U' 


3u' 


u. 


i  9t 


+  u .  u 


+  u.U' . 


i  j  9x' .  i  j  9x'  . 

J  J 


+  u. 


9u'  u'  . 
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and 


_  ^  n  !  ij  t 

el?l  +  eu’.  fT-  +  eu-.  |®L+  e— J- 


9t 


-  0 


j  9x'  . 
J 

30  fu! 
_ 1 

9x '  . 

J 


j  9x' . 
J 


9x'  . 


J 


2o  i 


v'0 


9^0 
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(3.18) 


AdcHri£^Eqs .  3.15  and  3.17,  and  Eqs.  3.16  and  3.18  respectively 
and  using  the  fact  that  the  turbulent  quantities  at  x  are  indepen¬ 
dent  of  the  coordinates — jl! — and  vice  versa  result  in 
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r.  =  x'.  -  x.  blST.  "BET.  (3.21) 

and 

xlij  =  J  (x'  +  x  )  ,  CENTROID  (3.22)  /.oCffTlav 


it  results  in 


3 

1 

3 

“t n 

3x . 

2 

3x"  . 

3r.  1 

3 

J 

_ 

(3.23) 
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(3.24) 

(3.25) 

(3.26) 


Using  the  above  transformations  3.21--3.26  and  taking  average 
of  Eqs.  3.19  and  3.20,  we  obtain 
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(3.27) 


and 
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In  order  to  make  Eqs.  3.27  and  3.28  simpler,  the  local 
homogeneity  i  <;  a  tinned.  for  the  turbulent  fields,  i.e.,  the^sjDatia^_ 
variation  o£the  averaged  turbulent  quantities  is  negligible  over  a 


distance  r  <<  £  ,  where  £  is  the  scale  of  energy-containing 

—  S  e 

eddies.  Thus, 
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make  Eqs.  3.27  and  3.28  become 
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in  which  the  averaged  turbulent  quantities  are  functions  of  r  and  t  . 

For  further  simplifying  Eqs.  3.30  and  3.31  without  loss  of 
generality  in  a  locally  homogeneous^turbulent  shear  flow  of  thermally 
stratified  fluid,  the  following  are  assumed 

v  =  v'  ,  VT  =  >  g3  =  g'3  =  g  >  g2  =  g'2  =  81  =  g'l  =J?  * 


IL  =  constant  ,  U_  =  constant  ,  p  =  p' 

2  3  oo 


U.  =  f(x_)  ,  T  =  g(x  )  , 


where  is  in  the  streamwise  direction,  is  in  the  lateral 

direction,  and  x7  is  in  the  vertical  direction. 


and 
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if  i  'Y  j  •  Expanding  UJ  into  a  Taylor  series  at  gives 


dU, 


'!  d2U  (X’  -x  ) 

ui  -  U1  =  d3^-  <X*3  -  X3)  +  ^4  ■ 


2! 


dll1 

SJ  r3 


(3.34) 


if  we  assume  that  r^  <<  (dU1/dx3)  /  (d2U1/dx23)  . 

Hence,  substituting  To,  3.34  into  Eqs .  3.32  and  3.33  and 
contracting  Eq.  3.32  gives 
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in  which  the  pressure  velocity  correlation  terms  are  eliminated  in 
the  contracted  tensor  form  because  of  the  incompressibility  of  flow, 
i.e.  , 


P'  ,u. 


9r .  1  i 

1 


P,u'.  =  0 


9r .  1  i 

l 


In  order  to  transform  Eas.  3.35  and  3.36  into  wave  number  k 
space,  the  following  definitions  are  given 


(3.37) 


(3.39) 


(3.40) 


(3.41) 


(3.42) 


(3.43) 


lim  E.  .  =  lim  E. 


k  ^°° 

1,1 

k  ’ 1 

V  k  -I 

k  ooT 

u_ 

3  • 

■L3 

3  • 

are  used.  It  is 


realized  that  E^  ^  and  E^  q  are  the  cospectra  of  the  Reynolds 
stress  and  vertical  heat  flux.  Since  the  above  two  equations  are 
expressed  in  k  space,  the  directional  information  of  the  spectra  is 
retained  in  both  equations;  needless  to  say,  it  is  difficult  to  solve 
Eqs .  3.42  and  3.43.  Now,  if  Eqs .  3.42  and  3.43  are  averaged  over  a  V 

spherical  shell  then  we  have  the  equations  which  are  only  function 
of  k  ,  the  magnitude  of  k  ,  and  the  directional  information  is 


J 


(3.44) 


(3.45) 


in  which 
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A 

1 


<Kk,t)  =  J  f  Ei4(k,t)do(k) 


♦w(k,t)  =  ^E1>3(k,t)do(k) 


3E  f  3E  (k,t) 

<kl  -Sk7  )Sp.  av.  1  2  i  kl  -^7 - 


"TT 


(k,t)  =  f  ET(k,t)do[k) 

A 


^WT(k  > t)  =  J  Mu  e  (k ,  t)dq  (k) 
'**  3 

F(k,t)  =  j  £  ^(lc.tJdoCk) 
FTT(k,t)  =  j>  F  (k,t)do(k) 


3E„ 


3E  (lc,t) 

(R1  air  )  =  p  ki  — ir — do(k) 

1  3k3  sp.  av.  3k3  - 


,  =  _J_  dU  _  dUl  dT  _  dT  _  dT 

—  9  dz  dx7  9  dz.  ^  dx7  dz 

l _ 1/  3  JU  3 

7 


(3.46) 


and  da(k)  is  the  surface  element  on  the  spherical  shell  with  radius 
k  . 


At  this  stage,  it  would  be  worth  describing  the  physical 
sense  of  each  term  in  the  spectral _Eqs_p  3«4fl  and  3.45  before  attempting 
to  solve  them.  The  first  terms  of  Eqs .  3.44  and  3.45  are  the  rate  of 
change  of  turbulent  energy  ajid  temperature  inhomogeneity  in  their 
corresponding  spectral  forms.  In  case  of  steady  turbulent  flow,  both 
terms  vanish.  The  second  terms  of  Eqs.  3.44  and  3.45  represent  the 
production  of  turbulent  energy  due  to  the  work  of  the  Reynolds  stress 
against  the  mean  shear  and  the  production  of  temperature  inhomogeneity 


35 


transferred  by  the  vertical  heat  flux  against  the  mean  temperature 


gradient.  The  integration  of  4>  and  <J>  with  respect  to  k  over 


the  internal  (0,  00 1  will  give  the  Reynolds  shear  stress^and  the 

vertical  h  e  a  t__ £  1  nx ly T 

The  third  terms  of  Eos.  3.44  and  3.45  are  the  turbulent  energy 
transfer  and  the  temperature  inhomogeneity  transfer  respectively  due 
to  distortion  bv  mean  shear.  Unlike  the  second  terms  of  Eqs.  3.44  and 
3.45,  these  terms  do  not  contribute  the  total  energy  and  temperature 
inhomogeneity  balances  but  redistribute  energy  and  temperature  inhomo¬ 
geneity  by  transfer  through  wave  numbers.  The  same  situation  happens 
to  the  fourth  terms  of  Eqs.  3.44  and  3.45,  i.e.,  these  are  also  the 
transfer  terms,  however,  these  transfer  terms  are  due  to  distortion  by 
the  fluctuation  gradients,  or  say,  due  to  ipertial  processes.  These 
inertial  transfers  from  low  wave  numbers  to  high  wave  numbers  are  cer¬ 
tainly  the  characteristics  of  turbulent  flow. 

As  to  the  last  terms  of  Eqs.  3.44  and  3.45,  they  are  energy 
dissipation  by  viscosity  and  temperature  inhomogeneity  smeared  out  by 
thermal  conductivity.  Now,  here  comes  the  most  important  term  j 

Jty^rr(k,t)  ^  in  Eq.  3^44  which  reflects  the  effects  on  the  turbulent 
spectra  by  buoyancy  force  due  to  stratifications.  In  the  case  of 
stable  stratification,  this  term  becomes  a  fsink  with  respect  to  the 
turbulent  energy,  and  on  the  other  hand,  it  becomes  an  energy  source 
in  case  of  unstably  stratified  flow. 

In  the  following .  _,the  turbulent  flow  is  assumed  to  be  steady 
state,  in  the  wave  number  range  investigated.  And  also  the  Ijhird  fterms 
in  Eqs.  3.44  and  3.45  are  assumed  to  be  negligible  as  compared  to  the 


fourth 


ones  because  it  is  believed  that  these  shear  transfer  terms 
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,yyy c 


mainly  affect  the  spectra ^iii^tkg^ower  wave  number  range  t h an_b uoy an c y 
subrange .  Thus,  the  spectral  equations  for  velocity  and  temperature 


fields  are  respectively  the  following: 


l: 


F(k)  ■  *uw(k)  HI  +  WwT(k)  "  2vk2^k^  =  0 


and 


whose  integrated  forms  are 


and 


(3.47) 


(3.48) 


(3.49) 


(3.50) 


In  Eqs .  3.49  and  3.50,  the  total  dissipation  of  turbulent 
energy  by  viscosity  and  N@  twice  of  the  total  dissipation  of  tem¬ 
perature  fluctuation  by  thermal  conductivity  are  defined  by 


c  =  2 v  /  k24)(k)dk 

o 


(3.51) 


and 


N*  =  2v?  /  k2  d>TT  Ck)  dk  =  _2N 
o 


(3.52) 


In  the  sections  following  the  generalized  eddy-viscosity. 

/  - — ■ — ■  *“*  1 

approximation  is  introduced  and  closed  forms  of  solutions  of  Eqs .  3.49, 
and  3.50  are  analytically  derived  for  different  situations. 
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3. 2  Generalized  Eddy-Viscosity  Approximation 

As  a  consequence  of  the  nonlinearity  of  the  Navier-Stokes 
equations,  the  correlation  equations  like  Eqs.  3.32  and  3.33  always 
contain  one  more  unknown  than  the  number  of  equations;  in  other  words, 
a  closure  problem  is  involved.  Thus,  if  the  problem  must  be  solved 
from  the  correlation  equations  or  its  corresponding  spectral  equations, 
additional  assumptions  have  to  be  assigned  to  the  turbulent  energy 
transfer  function  if  only  the,,  locally  isotropic  turbulence  is  con¬ 


sidered.  In  case  the  temperature  field  and  shear  flow  are  introduced, 

- - - - — 


i  otherjas sumptions  should  he  addp.dj  Thus,  the  six  unknowns  in 
Eqs.  3.49  and  3.50  are  reduced  to  two  unknowns  and,  of  course,  these 
equations  are  solvable. 

In  the  following,  the  generalized  eddy-viscosity  of  the  form 


(3.53) 


The  physical  sense  of  the  generalized  eddy-viscosity 
expressed  by  Eq.  3.53  is  not  difficult  to  be  realized  if  Heisenberg 1 s 
idea  (1948)  will  be  reviewed.  His  form  to  express  the  turbulent 
energy  transfer  function  is  listed  as 


/  F (k) dk  =  n^k)  /  2k2<Hk) 


dk 


(3.54) 


and  r)  (k)  called  the  kinematic  eddy-viscosity  is  expressed  as 

„  (I)  Ph 

n^k)  =  Yl  J  ^WkS-^dk 


(3.55) 


where  y ^ 


is  a  numerical  constant. 
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The  idea  implied  in  the  Eq.  3.54  is  that  the  energy  transfer 
from  the  wave  number  less  than  k  to  the  wave  number  larger  than  k 
can  be  considered  as  eddy-viscosity  working  on  the  turbulent  vorticity 
formed  in  the  wave  number  interval  (0,k).  This  eddy-viscosity  can 
be  viewed  as  the  integral  effect  of  eddies  with  wave  numbers  larger 
than  k  on  the  eddies  with  wave  numbers  less  than  k  .  Thus, 
according  to  the  dimen^  jr>na i  axgjmaaX g ,  n^(k)  is  expressed  in  a 

form  expressed  by  Eq.  Now,  we  can  see  that  the  expression  of 

the  generalized  eddy-viscosity  ng(k)  has  rp-rtaiplv  thp  dimpn^ipp  nf 
‘the  eddy-viscosity,  moreover,  the  parameter  s  introduced  in  Eq.  3.53 
can  be  interpreted  as  degrees  of-interaction  between  the  motions  of 
eddies.  i.e.t  how  the  motions  of  eddies  are  interrelated  to  one 
another.  Generally  speaking,  the  introduction  of  s  will  aot  affect 
the  spectral  form  but  its  magnitude  as  can  be  seen  later. 


As  will  be  seen  later,  in  case  s  ->  °°  ,  the  expression  of 

Eq.  53  becomes  only  a  function  of  the  local  wave  number  k  and  the 

associated  energy  spectrum.^  Of  course,  Heisenberg1  s  expression  for 

the  eddy  viscosity  is  only  a  special  case  of  Eq.  3.53  when  s  =  1  . 

^  *  —  - 

Based  on  the  same  argument  and  for  the  purpose  of  further  generaliza¬ 
tion,  the  turbulent  energy  transfer  function  is  expressed  as: 


/  F(k)dk  =  Y 

L.  ^  >  1 


s  s  .  . 

“  7  -  7  -  1 

/  <j>(k)k  dk 


, ~  r  3r-l 
/  2*  (k)k dk 


LO 


(3.56) 


which  can  evidently  be,  derived  from 


/  F(k)dk  =  const, 
k 


dl  d2 
(.  (k)k  dk 


~|d3 


k  d^  d 
<t>  (k)k  dk 


(3.57) 
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given  by  Goldstein  (1951) ,  where  the  exponents  are  related  from 

■  um  i  '■  ■— 


dimensional  arguments  as^ 


and 


d3(d2  *  1)  *  d6Cd5  *  1)  =  f  , 


Vl  *  d6d4  "  2  ' 


These  two  relations  can  reduce  six  parameters  d^  ^in 
Eq.  3.57  into  four  arbitrary  parameters;  for  further  simplification* 
Eq .  3.56  is  eventually  derived.  It  can  readily  be  seen  that  the 
\ First  factor  in  Eq.  3.56  has  the  dimension  of  kinematic  turbulent 


^  eddy  viscosity  and  the^second^  one  has  a  dimension  of  turbulent  vor- 
ticitv.  It  is  clear  that  the  expression  of  Heisenberg's  form  of 
turbulent  energy  transfer  Eq.  3.54  is  a  particular  form  of  Eq.  3.56 
in  case  s  =  1  and  r  =  1  .  _ 


S  i m i  larl\-,  the  transfer  of  turbulent  temperature  inhomogeneity 

00  * 

^  F^(k)dk  can  also  be  in  a  form^of^enerali^g^  eddv-viscosity *  i.e. 


- 

"  k  -  -  - 

/  2^TCk)k5r_1  dk 

r 

1 

-  0  , 

(3.58) 


where  b  is  a  numerical  constant  of  order  1  and  is  equivalent  to  the 
ratio  of  eddy  thermal  diffusivity  to  eddy  kinematic  viscosity. 

Before  expressing  (f^Ck)  and  4)uw (k)  in  terms  of  the 
generalized  eddy-viscosity,  the  validity  and  generality  of  Eq.  3 . 56 
will  be  seen  as  follows. 
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lim 

L  k 


Since 
s 

CD  - 

r  ,2 


(k)k 


-  =r  -  1 


dk 


1_ 

s 


J_  }_ 

2  2 

<t>  (k)k 


(3.59) 


and 


lim 

p-K» 


,r  3r-l 
t>  (k)k 


dk 


4>(k)k3 


Eq.  3.56,  in  case  s  ->  00  and  r  ->  00  ,  becomes 


f  F  (k)dk  =  y 

CD  ^  CD  °°, 


3  5 

2  2 

<t>  (k)kz 


(3.60) 


It  is  obvious  that  Eq.  3.60  is  nothing  but  Kovasznay's 
approximation,  which  is  exactly  the  local  limit  of  the  nonlocal  and 
generalized  apprnYimnt i  nn  n.f_  the,  fprm  of  Eg.  3.36.  From  the  example 
illustrated  above  to  indicate  the  usage  of  Eq.  3.56,  we  can  say  that 
by  varying  the  values  s  and  r  ,  solutions  corresponding  to  the 
different  degrees  of  turbulent  nonlinear  interactions  between  eddies 
can  be  obtained.  Obviously,  the  Kovasznay's  approximation  is  a  limit 
form  since  it  means  that  the  eddies  interact  themselves  only.  For  the 
case  0  <  s  <  00  and  0  <  r  <  <»  ,  it  may  be  interpreted  that  the 

eddies  of  wave  number  k  interact  with  the  other  eddies  with  wave 

number  k  ±  Ak  ;  of  course,  Ak  is  a  function  of  the  values  of  s 

and  r  ,  and  Ak  decreases  with  increasing  s  and  r 

Keeping  these  ideas  in  mind,  the  generalized  eddy-viscosity 

1  1  •  —  i  r 

approximation  can  be  extended  to  the  spectra  ^^(k)  and  4>wT(k) 

without  any  trouble.  In  the  following,  r  =  1  will  be  assigned 
r  rk  2 

since  [  J  4)(k)k  dk]  gives  a  clear  physical  sense--the  spherically 
averaged  square  of  the  root-mean-square  vorticity  of  turbulence_in 


41 


the  range  0  to  k  .  The  exact  local  approximation  when  s  -*  <»  and 


r  oo  will  also  be  used. 


1VJ  Jy* 

In  case  <J>  (k)  and  d>  ^(k)  are  concerned,  we  need  consider 

_  uw^  J  wT  _  —  ■  -  - - 

not  onl^the  turbulent  field  but  the  interactions  between  the  mean  , 

velocity  and  temperature  fields  and  the  turbulent  field. 

■—l — 


where  the  upper  sign  indicates  the  case  when  dU/dz  >  0  ,  and  the 

lower  sign  for  dU/dz  <  0  ,  and 


and  the  lower  sign  denotes  the  unstable  one  dT/dz*  <  0  because  the 
vertical  heat  flux  is  negative  in  case  of  stable  stratification  and 
is  positive  for  the  unstably  stratified  flow. 

The  parameters  in  Eqs .  3.61  and  3.62  are  reflated  from 

dimensional  argument  as  C,  +  CL  =  1  ,  and  C-  +  C.  =  1  .  C,  and 

- h —  i 1  j — /  v3  4  j  a - 

characterize  the  degrees  of  interaction  between  the  mean  velocity 
gradient  and  the  turbulent  vorticity;  and  denote  the  degrees 

of  interaction  between  the  mean  temperature  gradient  and  the  turbulent 
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temperature  gradient.  In  case  .C,  =  0  ,  C,  s  1.  ,  and  C.  =  1  ,  CL  =  0, 

- - - - -  '-J— - L — — *  Li _ i _ ? 

Eq.  3.61  expresses  the  conditions  of  strong  interaction  and  Weak  inter¬ 


action  of  velocity  field  considered  by  Tchen  [19531 .  Hence,  increasing  r 
s  would  mean  that  the  interaction  between  the  vorficitv  of  main 
motion  and  turbulent  vorticity  becomes  stronger,  and  according  to 
Tchen  (1953),  the  resonance  between  two  motions  is  intensified.  Thus, 
in  case  resonance  is  intensified,  the  inertial  transfer  process  through 
eddies  is  interfered  by  the  mean  motion  and  energy  is  supplied  to  the 
eddies  by  means  of  the  Reynolds  stress  working  against  the  mean  motion. 

In  other  words,  the  -  5/3  law  is  invalid  and  the  energy  spectrum  has 
,less  steep  slope.  In  the  extreme  case  when  C0  =  1  ,  the  slope  of 
the  energy  spectrum  becomes  -1. 

As  far  as  stratificat ion  is  concerned,  »  0  ,  =  1^ ,  and 

=  1  ,  =  0  ^are  equivalent  to  the  case  of  strpng  and  weak  inter¬ 

actions  investigated  by  Gisina  (19661.  In  the  following  sections, 
the  significance  of  will  be  seen.  Thus,  it  is  clear  that  Eqs.  3 . 6 1 

and  3.62_can  give  more  general  information  on  the structure  of  the  ther^ 


iMjjy  i  fied  turbulent  shear  flow,  since  different  values  of 

can  be  assigned  to  characterize  different  flow  conditj ons . 

T/  t(x.  *( 

In  ^ddition^to  the  relationship  between  C. 


cribed  pre¬ 
viously  it  seems  to  be  helpful  to  let  >  0  and  >  0  .  Since 

|dT/dz*|  0  would  mean  that  the  flow  becomes  lesser  stratification, 
the  flow  will  be  nonstratified  when  |dT/dz*|  =  0  ,  thus,  it  implies 

that  there  exists  no  vertical  heat  flux.  But  if  =  0  ,  this 

implication  cannot  be  seen  when  |dT/dz*|  =0  ;  in  other  words,  if 

=  0  is  assigned  to  Eq.  3.62,  we  always  have  the  vertical  heart, 
flux  even  when  |dT/dz*|  =  0  .  However,  once  |dT/dz*|  ^0  ,  from 
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the  dimensional  arguments  can  be  less  than  or  equal  to  0 

according  to  the  relation  =  1  .  Thus,  it  may  conclude  that 

is  a  function  of  |dT/dz*|  ,  i.e.,  in  case  |dT/dz*|  -*•  0  ,  it  is 

L  £  i  —  —  ■  ■ 

believed  that  C7  must  be  greater  than  0  as  stated.  Similar  argu- 
--  ^  — 1 1  " 


ment 


can  be  applied  to  if  Eq.  3.61  is  considered. 


3^3  Solutiflps  of  the  Spe^ji^al  Equations  of  Thermally  Stratified 
Turbulent  Shear  Flows 

In  this  section,  the  analytical  derivations  of  solutions 
based  on  the  generalized  eddy-viscosity  approximation  are  described 
in  detail.  Nondimensionalized  spectral  equations  are  derived  not 
only  because  they  can  give  neat  and  concise  forms  but  because  some 
similarity  theory  of  spectra  can  be  made  if  suitable  dimensionless 
variables  are  used. 


3.3.1  Solutions  of  flows  wrfch,  molecular  effects  -  The 
introduction  of  molecular  effects  does  not  mean  that  the  spectra  at 
high  wave  numbers  can  be  studied  by  the  present  method  as  seen  from 
the  discussion  in  section  2.1.  Thus,  only  the  spectra  in  the  buoxan&y# 
and  inertial  subrange  are  of  interest. 

Based  on  the  generalized  eddy  viscosity  approximation 
described  in  section  3.2,  Eqs.  3. 49  and  3_, 50  together  with  Eqs.  3.56, 
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45 


energy  and  temperature  spepfTf1  ('an  ha  aiLalualAcLky 
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(3.72) 


(3.73) 


As  to  the  turbulent  shear  stress  spectrum,  and  the  vertical 
heat  flux  spectrum,  Eqs .  3.§X  and  3.62_  are  transformed  into  non- 


dimensional  forms.  Using  the  dimensionless  variables 


r 


/$  (x)dx  =  - 
J  uw  J 

X 


1  -  Y“ 


|  m  |  Y  +  Y  2  +  g  Z  4 
1  1  s  s  Is 


and 
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/  <I>wT(x)dx  =  + 


1  -  Y 


m 


|Y  +  Y  2  +  6.Z  4 

1  s  s  Is 


Ca  s 


fications,  the  spectral  equations  differ  only  in  the  turbulent  energy 
spectral  equation  as  indicated  in  Eq.  3.68.  Thus  the  turbulent  energy 
spectrum  is  explicitly  influenced  by  thermal  stratification  while  the 
temperature  spectrum  is  implicitly  affected  due  to  the  introduction 
of  energy  spectrum  in  a  form  of  eddy  viscosity  as  shown  in  Eq.  3.69. 

Using  Eqs.  3.6J5  and  3.69^,  a  relationship  between  Ys  and  Zg 


can  be  deduced  from 


C. 

I m I Y  2  +  Y  2  +  g.Z  4 
1  1  s  s  1 


Pr<zs2  *  MV) 


(3.79) 
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Thus,  Eq.  3.79  together  with  either  Eq.  3.68  or  Eq.  3.69,  and  Zs 

can  be  solved  explicitly  or  implicitly  as  functions  of  x_ ,  and  con¬ 
sequently,  4>(x)  ,  <I>^,(x)  ,  4>uw(x)  ,  $w^,(x)  can  be  evaluated 
numerically  by  using  Eqs.  3^7^>  3.73,  3.77  and  3 . 7 8_.  However,  Tn 
case  C^= 1  ,  C^=0  ,  some  analytical  closed  forms  can  be  obtained. 

For  the  sake  of  interest,  the  following  analytical  solution  is  derived. 
First,  from  Eq.  3.79  and  ^^=1  ,  we  have 


C  _  C 

[(Pr-l)Zs2  +  |mT|PrZsa  ♦  l]Ys2  +  |m|(l-Zs2)Ys  +  (1-Zs2) 


Pr(Zs2  +  |mT| Zs4)  =  0 


(3.80) 


Now  let 


A(Z  )  =  (Pr-l)Zs2  +  ImJPrZ  4  +  1 


(3.81) 


B(Zs)  =  |m|(l-Zs2) 


C  C 

C(ZS)  =  +  BxZs4  (1-Zs2)  -  Pr(Zs2  +  |mT|Zs4)  , 


(3.82) 

(3.83) 


then , 


-B  + 


A/B2- 


4  AC 


2A 


(3.84) 


..  .  v  -B  -  V  B2-4AC  .  .  | „  i 

the  other  root  Y  =  - Z- - is  omitted  since  Y 

s  2A  1  s 1 

less  than  or  equal  to  1  and  Y^  ->  0  ,  Zs  -+  0  as  x  ->  0 


must  be 


Squaring  Eq.  3.84  yields 


Y  2 
s 


2B2-4AC-2B ~V  B2-4AC 


4A2 


(3.85) 
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Differentiating  Eq.  3.85  with  respect  to  x  gives 


,  dY  2 

Y  Y '  =  —  - 

s  s  2  dx 


1 

8A4 


A2  {  4B  41-  Z '  -  4A  4~  Z*  -  4C  44“  Z^  -2  44-  Z^  V  B2-4AC 


dc 


dA 


dB 


dZ  s 
s 


dZ  s 
s 


dZ  s  dZ  s 
s  s 


V 


- -  (2B  4|-  Z'  -4A  4f-Z'  -  4C  4^~  Z') 

-  v  dZ  s  dZ  s  dZ  s 

s  s  s 


B2  -4  AC 


-  2  j  2B2-4AC  -  2B  V  B2-4AC  j  A  Z ' 

s 


S  )A2 


8A4 


b§_.a£-.c^ 

dZ  A  dZ  dZ 

S  S  S  / 


2~\f  B2-4AC  4| - r-B— '■ 

dZs  Vb2-4AC 


dB 


dC  dA 


28  dT-  "  4A  IT  -  4C 


dZ 


s  J 


-2  J  2B2-4AC  -  2bV  b2-4AC  )  A 


dA 

dZ 


in  which 


dA  i  i-  ^4”^ 

ar  =  2(pr-1)zs  *  lmTlpr  c4  zs 

s 


(3.86) 


dB 

dZ 


=  -  2  m  Z 


and 


§-  -  *  Ws*  ‘  .  Bl(C4t2)z'4n-2PrZs-Pr|n,T|C4z'4  ' 
s 
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Also  from  differentiating  Eq.  3.69  with  respect  x  ,  we  have 


1-Z  2 
s 

s-‘  C^l 

1^1 

ZS4)("2ZS)' 

-(1-Z|)(2Zsh.C4| 

1^1 

Pr(Z2+ 

lraTl 

lzs^ 

Pr (Z2+ 

l"TlZs4)2 

£  3s  . 

?  i  ~  ^ 

=  -  (Y  Y '  j  x 
s  s' 

Thus,  from  Eq.  3.86 

„i>s  [|»vl(2-c4)Zs  *  2Zs  *  Cjm^4  ’] 

slprJ  c 

(Zs2  *  Klz/)5*1 


(3.87) 


4  B 


17  .  is 

|._  dB  ..  dC  dA  \ 

2B  dz  -4A  dz  -  4C  dZ 

S  s  s  .. 

-  2 

2B2-4AC-2B V  B2-4AC 

A^  1 
dZs' 

2 

(3.88) 
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Z'  =  x 
s 


3s-2  2  2s  2s-2 

2-5  ss_2(Pr)2"s(i-zs2)  2-5  rz  2  + 


2s+2 


(Zs2  *  KIZs4)  2"s 


C.  +  i  C  -1-1 

|mT|(2_C4)Zs  +  2Zs  +  C4|mT,Zs 


s-2 

I  1  V 

< 

00 

s 

2-s 


4  (B  S5--  A  -  C  %■ 
\  dZ  dZ  dZ 


-  2  V  B2-4AC 


dB  B 

dZ 


s  V  B2-4AC 


/  D  dB  dC  Ar  dA  ] 

(2B  dz  -  4A  dz  •  4C  dz  ) 

-  2 

2B2-4AC  -  2B V  B2-4AC 

A§-] 

S  S  s 

Si 

2-s 


(3.89) 

In  case  s=2  ,  the  abridged  equation  can  be  obtained  from 


Eq.  3.88, 


4'  3 


x4  = 


Pr2(Zs2  +  |mT|Zs^) 


C  .  + 1  C  -1 

16Ani-Zs2)[|mT|(2-C4)Zs4  +  2Zs  +  C4|mT|Zs4  ] 


4|B  %r  - A  §-  - c  -  2Vb2-4ac  f- 

s  s  s  s 


V  B2-4AC  '  dz 


2B  4^ - 4A  4S - 4C  dA 


dZ  dZ 

S  S  S  J 


2B2-4AC-2B V  B2-4AC 


dA 


dZ 


(3.90) 

and  in  case  2  <  s  <  °°  .  a  closed  integral  form  can  be  derived  from 


Eq.  3.89, 
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4s 

s-2  s-2 

—7—  x 

4s 


2  2s 

s2_S  (Pr)5"2  J (1-Zs2) 


2s-2 


s-2 


(2-C4)Z 


C4+1 

s 


+  2Z 

s 


+  C 


.Ki 


,C4-1 


s-2 


(Z, 


|i\l 


Zs^ 


2s+2 

s-2 


1 

8A4 


s 

i  s-2 


B 


dB 

dZ 

s 


2  V  B2-4AC  g- 
s 


B 

Vb2-4AC 


on  dB  ..  dC 

2B  Hz"  -  4A  dZ“ 
s  s 


4C 


dA 

dZ 


2B2-4AC-2BVB2-4AC 


dZ 


S 

s-2 


dZ  +D. 
s  1 


(3.91) 


where  the  constant  can  be  evaluated  from  the  boundary  condition 

x  =  0  ,  Z^  =  0  ,  or  x  =  °o>Zs  =  l  ;  the  integral  is  an  indefinite 
one.  Therefore  from  Eqs.  3.72,  3.73,  3.77,  and  3.78,  spectra  can  be 
obtained  as  functions  of  s  ,  Pr  ,  m^  ,  ,  m  ,  and  6^ 

3.3.2  Solutions  of  flows  with  negligible  molecular  effects  - 
In  the  previous  statement,  molecular  effects  as  represented  by  the 
kinematic  viscosity  and  the  thermal  diffusivity  are  involved.  How¬ 
ever,  in  case  the  spectra  at  waves  numbers  far  away  from  the 

Kolmogorov's  wave  number  is  considered,  i.e.,  when  k  <<  ,  the 

^ 

local  dissipation  _of^_ turbulent  energy  and  the  local  dissipation  of 
temperature  inhomogeneity  in  the  range  0  to  k  are  negligible  as 

compared  to  and  N  ,  respectively.  This  will  be  the  case  when 
we  consider  the  buoyancy  subrange  where  the  molecular  effects  are 
negligible.  Thus,  Eqs.  3.49  and  3.50  are  reduced  to 
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(3.92) 


(3.93) 


Again,  the  generalized  eddv-visnnslt.v  annrnyimat-i  rm  is  applied; 
thus,  we  have 


(3.95) 
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Now,  with  the  nondimensional  variables 


(3.96) 


(3.97) 


(3.98) 


and 
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Thus,  from  Eqs.  3.98,  3.99,  3.101,  and  3.102,  4>(x)  ,  4>  (x)  , 

n— — — W  TT 

4>  (x)  ,  and  <J>wT(x)  can  be  evaluated  numerically.  However,  for 

simplifying  the  analytical  derivation  of  solutions,  two  sets  of  values 


of  s  and  r  are  assigned. 
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(A) 


2  <  s  <  00  ,  r  =  1 


By  letting 


x  Y  (x)Y'(x) 

Y  2(x)  =  /  2x2$(x)dx  ,  or  4>(x)  =  — - - - 


(3.103) 


and 


x  Z  (x)Z'(x) 

Zc2(x)  =  /  2x2$TT(x)dx  ,  or  4>TT(x)  =  — - — - 


TT 


(3.104) 


Bqs.  3.98  and  3.99  become 


/  (Y  Y')2  x  2  dx 

J  S  S 


*iT 

1  ^  A  i 

y  2  +|r  y  2  +  r  z  4 

s  1  1  s  Is 


and 


1  = 


s  3s  .  - 

OO  -  _  -  _  \ 

I  (Y  Y')2  x  2  dx 

J  s  s 

X 


(3.105) 


(3.106) 


Also  from  Eqs .  3.101  and  3.102,  we  have 

(V  *|r|ys2  *  qZ;4)  C2  Y^2  1  Y' 


4>  (x)  = 

uw 


(Y  2  +  | T  | Y  2  +  r.z4)2 

v  s  1  1  s  1  s' 


C  C  -1  C  -1 

Y  2 (2Y  Y»  +  I r I C0Y  2  Y»  7  r,C  7  4  Z'  ) 
svss  2s  S  14s  S 


C  C 

(Y  2  +  | r | y  2  +  r.z  4)2 

v  s  1  1  s  Is 


(3.107) 


and 
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\,T(X) 


C  C  -1  £  C  -1 

(zs2+|rT|zs4)c4zs4  z^-zs4(2zsz^  +|rT|c4zs4  z^) 

<zs2*  lrrlzs4’2 


(3.108) 


Now,  from  Eqs.  3.105  and  3.106,  we  have 


since,  by  differentiating  Eq.  3.105  with  respect  to  x  ,  x  can  be 

in  terms  of  Y  ,  Y'  ,  Z  ,  and  Z'  and  Eq.  3.109  can  offer  the 
s  s  s  s 

relationship  among  Yg  ,  Y^  ,  Zg  ,  and  Z^  .  However,  if  some  analyti- 


can  be  assumed.  Hence,  Eq.  3.109  becomes 

C 


y  2  +  r  y  -  z  2  -  (  rT  ±  r.)z 

s  1  1  s  s  1  T‘  Is 


4  =  0 


(3.110) 


and  so 


Y  = 
s 


- 1 r |  +  \/r2  +  4zs2  +  4 ( | rT |  ±  r1)Zg 


(3.111) 


Squaring  Eq.  3.111  results  in 

C. 

2r2  +  4Z  2  +  4(|r„|  ±  rjz  4-2|r| 

2  _  _ s _  1  T  ' _ I  s  11 


r2+4z2  +  4(|rT|±r1)zs 


(3.112) 


Inserting  Eq.  3.112  into  Eq.  3.106  and  differentiating  the 
resultant  equation  with  respect  to  x  ,  we  have  thus 
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an  abridged  form  can  be  seen  from  Eq.  3.113,  i.e., 


59 


4x4  =  (Z 


s 


2 


+ 


c 


lrTlZs 


3 


t2zs  +  c4lrrlV  1 


-1 

) 


C  -1 

2Zs+C4(|rTl±ri)Zs  '  |F 


2Vc4tlrTl±ri)s  * 


r2+4z2+4(|rT|±r1)Z5 


(3. 116) 


In  this  case,  an  exact  local  limit  of  the  generalized  eddy- 
viscosity  approximation  can  be  formed  and  Eqs .  3.98  and  3.99  become 


which  are  the  simultaneous  nonlinear  equations  of  $(x)  and  (x) 

  IT 


at  a  given  dimensionless  wave  number  x  .  Clearly  the  numerical  solu- 
tions  can  be  obtained  as  a  function  of  x  and  the  parameters  ^|r|  , 
Ti  ,  |  rT  |  ,  and  .  Since  both  parameters  C0  and  are 

retained  in  Eqs.  3.117  and  3.118,  and  the  spectral  equations  are  pre¬ 
sented  in  a  clear  and  simpler  form,  Eqs.  3.117  and  3.118  will  be 
investigated  intensively.  As  to  evaluating  $^(x)  and  $^(x)  , 
Eqs.  3.101  and  3.102  will  be  used.  Thus,  ^  ^ 
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(3.119) 


(3.120) 


3.4  Asymptotic  Solutions 

In  Previous  sectipns^_detail£ ixave to  solving 

the  spectral  equations  of  thermally  stratified  tnrhnlpnt._ shpar  flows: 

of  course .  numerical  s o lut i ons  are  avail ab le  now .  However,  if  we 

need  examine  the  significances  of  all  the  parameters  such  as  s  ,  , 

C2  ,  Cj  ,  and  introduced  in  our  generalized  eddy-viscosity 

— 

approximation,  it  would  be  helpful  to  investigate  the  as^TJHtotic 
solutions  under  certain  conditions.  In  order  to  avoid  the  tremendous 
complexity  introduced  by  the  consideration  of  molecular  effects,  we 
still  prefer  to  consider  the  case  when  the  molecular  effects  are 
negligible.  Evidently  in  the  buoyancy  subrange  of  wave  numbers  this 


is  the  case. 
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3.4.1  Asymptotic  solutions  of  stably  stratified  flow  -  First, 
consider  the  dimensionless  wave  number  range  x  >>  1  or  equivalently, 


1/2  3/4  v  3/4  -5/ 4  _3/2 

k  >>  y  b  N.  c  6 

s  ,r  * 


c  .  In  this  wave  number  range, 
the  buoyancy  effect  on  the  motions  of  eddies  is  negligible  as  com¬ 


pared  to  the  inertial  interaction  among  eddies  as  we  can  prove 
aposteriori  in  the  next  section.  Thus,  it  will  not  be  surnr-isina  that 
the  classical  -5/3  law  holds  in  the  dimensionless  wave  number  ran^e 
x  >>  7”  From^Tq! 3  1 T5^)  asymptotic  solutions  are  available  by 
assuming  |  F  j  "  0  ,  and  Z^2  >>  |r\p|Z  ^  as  follows: 


s-2  4-s 


,  3s  -  3s  3s  3  ,  0  -i 

3  2  s  x  A  2  <  s  <  «  ,  x  »  1  , 


(3.121) 


(3.122) 


Thus,  equivalently,  the  corresponding  velocity  and  temperature 
spectra  are 


(3.123) 


As  to  s  =  2  ,  Lq.  3.116  or  the  limit  form  of  Eqs .  3.123  and 
3,124  as  s  2  can  be  used,  and  we  have  then  ( Mi  ) 


X 


(3. 125) 
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s  =  2 


x  >>  1 


(3. 126) 


In  case  x  <<  1  ,  the  situation  is  somehow  more  complicated 
since  the  buoyancy  effects  will  distort  the  ijTertial  interaction  among 
eddies.  Two  cases  will  be  considered  in  the  dimensionless  wave  number' 
range  x  <<  1  .  In  both  cases,  we  still  assume  that  the  local  pro¬ 

duction  of  turbulent  energy  in  the  wave  number  range  considered  is 
negligible .  From  the  computational  point  of  view,  we  can  let 
in  Eq.  3.115.  Thus,  with  the  further  assumption 
lrTlzs4  >:>  z/  >  Ecl-  3-115  gives: 


I!) 


s-2 

3s 


,  |  _  |  _  ,  3  3s 
( | r? I +r1)  s  x 


2  <_  s  <  00 


and  from  Eq.  3.109, 


Y2  = 
s 


2  <  s  <  »  .(3.127) 


Thus ,  the  velocity  and_ temperature  spectra  become  in_this  case 
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These  asymptotic  solutions  can  be  considered  as  the  spectra  in  the 
inertial  and  convective  subrange  where  the  in^^^l  transfer  process 
predominates  the  velocity  field  and  the  interaction  between  the  main 
and  turbulent  temperature  fields  plays  the  principal  role  in  develop¬ 
ment  of  temperature  fluctuations. 

Next,  some  asymptotic  solutions  in  the  buoyancy  subrange  will 
be  investigated.  In  the  buoyancy  subrange,  the  local  production  of 
turbulent  energy  and  turbulent  temperature  inhomogeneity  is  so  small 
that  the  inertial  transfer  processes  predominate  turbulent  motions. 
However,  the  buoyancy  subrange  differs  from  the  inertial  subrange  in 
the  fact  that  in  the  buoyancy  subrange  the  buo^anc^  effects  due  to 
the  vertical  heat__flux  affect  the  energy  balance.  Thus,  with  the 
assumptions 

and  X 


1  r  1 

~  o  , 

^4  9 

r.z  4  >>  z2  >> 

|  r  |  z  4 

r,  >> 

\LL 

— / 

Is  S 

1  T1  s 

l 

C  -1  C  -1 

c4(lrTl  +  rpzs4  »  2zs  »  c4irT|zs4 


Eq.  3.115^gives  the  asymptotic  solution  as  follows: 


(3. 130) 


<<  x  <<  1  , 


C.  >  0  and  »  >  s  >  2 
4  — 


(3.131) 
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or  equivalently  we  have  the  asymptotic  temperature  and  velocity 


C  >  0 
4 


and  °°  >  s  >  2 


(3. 133) 


Now  the  significance  of  can  be  seen  clearly  from  Eqs , 


and 


3.132.  and  3.13£.  In  case  |C^  ->•  0  and  lj  ,  $(x)  ~  x' 

<I>  (x)  ~  x  ^  for  x  <<  x  <<  1  ,  and  this  is  exactly  the  solution 

H-— -  e 

predicted  by  Lumle^£1964} .  If  the  expression  (3.62)  for  the  vertical 
heat  flux  spectrum  is  reviewed,  0  and  -*■  1  would  mean  f±Lat 


the  vertical  heat  flux  spectrum  is  determined  bv  thp  velocity  field. 

oo 

implied  by  the  eddy-viscosity  y  ^  $S^(k)k~  s^‘'"  ^dlc]  in 

Thus,  after  reviewing 


Eq.  3.62  and  is  proportional  to  |dT/dz*| 

Lumley-Shur ' s  hypothesis  in  section  2.2.3,  the  generalized  eddy- 
viscosity  approximation  considered  at^  -»  0^  ,  ->■  j,  ,  [r|  ~  0^ 

Tl  is  equivalent  to  Lumley-Shur '  s  l|vpothes.i  s  .  However, 


and  T,  >>  T, 


1 


as  we  can  see  later,  because  of  the  generality  of  the  generalized 
eddy-viscosity  approximation  implied  by  varying  ,  C0  ,  , 


r|  and  r.  ,  the  present  results  are  more  fruitful.  And 
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5, 135,  then  ^>(x)  ~  x~  and  |$TT(x 


)  ~  x  j  are  obtained. 


This  is  actually  the  case  considered  by  Monin  (1962)  and  Gisina  (1966) 


3,4.2  Asymptotic  solutions  of  unstably  stratified  flow  -  In 
case  the  unstable  stratification  is  concerned,  it  is  better  to  study 
Eqs .  3.117  and  3.118  instead  of  Eq.  3.113.  For  the  same  reason  des¬ 
cribed  above,  let 
Thus ,  we  have 


r  I  ~  0  and  >>  |  T  I  in  the  buoyancy  subrange. 


3f  f 

3  3  1  ^4  1 

2  2  ”2  2  2  2 

1  =  x  0  (x)  +  TjX  *  (x)  $TTZ  (x) 


(3.126) 


and 


1  =  x 


$  (x)$TT(x) 


(3.127) 


In  the  wave  number  range  x  >>  1  ,  the  inertial  transfer  pro¬ 
cess  dominates  the  turbulent  structure T  it  can  be  expected  that  the 
-  5/3  law  hold$  for  both  velocity  and  temperature  spectra.  However, 
where  x  <<  1  and  is  still  far  away  from  the  range  ^x^  in  which  the_ 
production  of  the  turbulent  energy  and  ^t^perature^JjibomQ^eneitv  pre¬ 
dominates ,  we  have  the  buoyancy  subrange  for  an  -unstably  stratified 
turbulent  flow,  and  the  asymptotic  solutions  would  be 


jfrCxT 


and 


for 


x  <<  x  <<  1 
e 


(3.128) 


)  /; 


,W 


This  shows  that,  in  an  unstably  stratified  turbulent  flow,  the  velocity 
field  absorbs  energy  converted  from  the  temperature  field  in  the 


O 
A J 

s 

r 

f\ 

u 

F 
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buoyancy  subrange,  and  presents  a  hump  in  the  velocity  spectrum.  This 
is  actually  what  we  expected  in  our  search  for  the  basic  mechanism  of 
unstably  stratified  turbulent  flow  in  section  2.2.1. 


3.5  On  Bplgiatyo’s  and  Liimley-^fi^y 's  Hypotheses  and  the  Modified 
Hypotheses 

After  introducing  asymptotic  solutions  by  the  generalized 
eddy-viscosity  approximation  as  presented  in  section  3.4,  it  would 
be  worth  comparing  Bolgiano's  and  Lumley-Shur 's  hypotheses  for  stably 
stratified  flow.  Phillips  (1965)  made  comparisons  between  these  two 
hypotheses.  In  his  paper,  Phil lips  classified  these  hypotheses  as 
follows : 


/  Lumley-§hur ? s  hvuothes 

1 Jl  "The  statistical  properties  of  the  components  of  the  turbu¬ 
lence  with  wave  number  k  in  the  inertia-buoyancy  subrange,  including 
the  energy  spectrum  and  the  buoyancy  flux  spectrum,  are  determined  by 
the  spectral  kinetic  energy  flux  e(k)  (^T)this  wave  number  and  not 


at  distant  wave  number  in  either  direction," 


2/  "The  spectrum  of  the  buoyancy  flux  in  ^hy^&ical  space  of 


stably  stratified  environment  is  proportional  to  the  mean  buoyancy^ 

*7  j  - 

gradient."  ^ 


/V 


^Bolgiano 1  s  hypothesis---*^ 

1.  In  the  buoyancy  subrange,  the  statistical  properties  of  the 
motion  are  determined  by  N  ,  g/T  and  wave  number  k  alone. 

Phillips  agreed  with  Lumley  on  Lumley-Shur ' s  second  hypothesis , 
tjut  did—flot  agree  with  the  first  since  it  is  still  doubtful  if  the 


Kolmogorov’s  hypothesis  can  be  extended  to  the  buoyancy  subrange.  As 
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to  Bolgiano's  hypothesis,  Phillips  stated  that  the  quantity  N  is 
certainly  a  property  of  the  turbulence;  it  is  a  local  property  of  the 
inertial  subrange  and  an  integrated  property  of  the  buoyancy  and  energy- 
containing  ranges,  but  it  is  not  a  local  property  of  the  buoyancy  sub- 
range  itself  as  can  be  seen  from  the  equation 


e(k) 


/  £L  '-1 

dz 
_g_ 

\  j 


N(k)  =  e  + 


dT  rl 
dz 


-S- 

\  j  / 


(3.134) 


which  is  derived  according  to  the  assumption  that  the  production  of 

turbulent  energy  is  less  important  in  the  buoyancy  and  inertial  sub- 
ranges.  (Phillips,  1965). 

Certainly  Phillips"  argument  is  right  for  a  general  stratified 
flow.  However,  as  stated  in  sections  2.2  and  2.3,  Lumley's  solution 
would  be  meaningful  only  when 


M 


1  \ 
dz  1 

2 

<<  1 

gN 

\  i 

(2.19) 


Gisina  obtained  the  same  restriction  fpp  the  existence  of  Bolgiano's  4- 

.iS, 


solution  as  expressed  in  Eq.  2, 


Also  in  the  present  study,  asymp- 


r|  -  0  and  F^  >>  I r  |  as  described  in  section  3.4. 


totic  solutions  for  the  buoyancy  subrange  are  obtained  with  the 
assumption 

After  we  review  the  definitions  of  F^  and  as  expressed  in 

Eq.  3.96,  r  >>  | rT|  is  certainly  equivalent  to 

i  dT  , 


'dz. 


*  _ 

T 


for  N*  =  2N 


-  <<  1 

•  °r  ;  'd;i  ~  i 

HH 

(3.135) 

T 

dT 

-  2  dT 

’  dz* 

"  2  dl  * 
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(  v  1/2  — 

which  means  that  the  internal  frequency  (N/c)  g/T  of  turbulence 

of  stably  stratified  flow  must  be  greater  than  Brunt-Vaisal a  frequency 

—  1/2  1/2 

of  the  external  flow  (g/T)  (dT/dz)  and  which  is  exactly  the 

same  expression  as  Lq.  2.19  derived  by  Lumley  and  Eq.  2.20  obtained 
by  Gisina. 

By  combining  Eqs .  3.1^  and  3 . 135 .  we  will  see  what  Eq.  3.135 
will  really  ^mply  in  the  buoyancy  subrange.  Thus,  rewriting  Eq.  3.134 
as 


dT 


dT 


Jz  e(k)  N(k)  _  £  dz 

Ng  e  N  Ng 


(3.136) 


and  assuming  e(k)/e  to  be  finite,  Eqs.  3.135  and  3.136  will  imply 
N(k)  z  N  in  the  buoyancy  subrange.  In  other  words,  the  objection 
I  from  Phillips  to  Bolgianofs  hypothesis  can  be  released  if  the  in- 
equality  of  Eq.  2.19  or  Eq.  3.135  can  be  accepted  as  the  basic 


/ 


4 


requirement  for  the  existence  of  buovancv  subrange.  As  we  can  see 
later,  irftcase  3.135  docs  not  exist  for  a  certain  flow,  .condition , 
As  or  equivalently,  the  condition  >>  1  1^  cannot  be  fulfilled,  there 

exists  no  buoy anc/_jjubr ang e^  Thus,  one  word  can  be  added  that  if 


there  exists  any  buoyancy  subrange.  N(k) _ ~  N  is  still  a  local  prop- 

erty.  Of  course,  in  the  energy-containing  range,  this  argument  breaks 


» erty. 


down  without  anv  doubt. 

Before  the  author  elaborates  on  the  modified  hypothesis 
proposed  for  the  stably  stratified  turbulent  flow,  it  is  worthwhile 
to  review  Lumley* s  idea  concerning  the  production  of  turbulent  energy. 
In  his  paper,  Lumley  ( 1965)  considered  the  effect  of  the  turbulent 
energy  production  on  the  turbulence  spectra  by  expressing  the  shear 
stress  spectrum  as 
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(3.137) 


(3.138) 


where  R^.  is  the  spectral^fluJ^^ichardson^jiuinbe^ . 

As  stated  by  Lumley,  in  Eq.  3.138,  both  production  and 
buoyancy  spectra  have  the  same  form;  thus,  there  is  no  range  of  wave 


4 


number  in  which  production  ^.s  unimuflyf  apfr .  ‘  However j  i£  Eqs  .  3.61  and 
3.62^are  used  $nd)  if  the  corresponding  values  ,  C0  ,  ,  and 

are  known  for  a  given  flow  condition  such  as  stratification,  mean 
velocity  gradient,  etc.,  the  difficulties  induced  by  Lumley^  model 
can  be  removed.  That  is,  by  varying /'(T.  ^  there  may  exist  a  range  of 
wave  number  in  which  unjjn^orta^  and  the  buoyancy  force 

certainly  ^dominates  the  flow.  Now,  if  Eqs.  3.137  and  3.61  are  com¬ 
pared,  the  former  corresponds  to  the  latter  when  -*  1  ,  C0  0 

i.e.,  when  the  shear  stress  spectrum  is  proportional  to  dU/dz  the 
mean  strain  rate.  This  condition  can  only  be  fulfilled  when  the  mean 
strain  rate  is  small  compared  to  the  eddy  strain  rate  (Tchen,  1953). 
Thus ,  the  validity  and  the  generality  of  Eqs.  3.61  and  3.62  certainly 
offer  a  better  opportunity  for  the  study  of  buoyancy  subrange  of  a, 

thermally  stratified  flow. 

Consequently,  the  (modified  hypotheses  If  or  a  stably  stratified 
turbulent  flow  are  given  as : 

( l y/  The  buoyancy  subrange  of  a  stably  stratified  turbulent 
flow  exists  when  the  local  production  and  local  dissipation  of 


turbulent  energy  is  not  important  in  this  wave  number  range  and  when 
the  flow  conditions  satisfy  the  criterion 


l-l 

dz 


N 

T 


<<  1  ,  or 


J_  ldTl 
Y  'dz1 


N 

e 


JU 

T 


(3.139) 


y  In  the  buoyancy  subrange,  the  spectra  of  turbulent  energy 
and  temperature  fluctuation  are  determined  by  N  ,  e  ,  g/T  ,  dT/dz  , 


<f>  (k)  =  N 


2C  +4 

4 

-4  +  4C  . 

4 

8-4C, 

4 

4+C  . 

4+C„ 

\  4+C  - 

4 

4 

e 

g 

4 

T 

/ 

C  - 12 
4 

4+C  . 


f1(s,rT,r1,c4)k 


2+2C  .  2-2C  . 

4  4 

4+C  4+C 

<j)TT(k)  =  N  e 


-4+2C , 


t  4+CL 


-3C.-4 

4 

4+C  . 


T  i 


f2(s,rT,r1,c4)k 


(3.140) 


(3.141) 


L 


where 


"1+C4  ^ 

=  N*  e 


"1+C4  -  UC4 

■  dT  | 

'  dz  ' 


rT  -  n^4  e2 *-c^u 

^ Li 


-2+C .  1-C. 

4  i  dT  i  4  dT 


dzJ 


dzJ 


C  >  0  °°  >  s  >  2  N  =  2N  LL_  =  ?  — 

L4  ’  b   £  »  *  z  ’  dz .  dz  ’ 


and  f ^  and  f2  are  numerical  variables  as  function  of  s  ,  T_ 


r,  and  C  . 

1  4 


In  case  C4  =  1  ,  thus  =  1  ,  then  Eqs .  3.133  and  3.134 

become  the  Bolgiano's  solutions  FQg  2  4  arlH  2  1 
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Now,  from  Eq.  3.140,  we  can  see  that  the  reason  why  the  parameter  e 
is  dropped  from  the  Bolgiano?s  hypothesis  is  not  because  e  <<  e(k) 


as  interpreted  by  Bolgiano  ( 1959)  but  because  =  1  makes  zero 
exponents  for  e  in  Eqs.  3. 140  and  3.141.  Accordingly,  the  parameter 
e  must  be  retained  in  the  hypothesis  as  described  in  the  hypothesis 
(2).  As  to  the  physical  background  on  which  the  modified  hypotheses 
^are  based,  the  section  4.4  in  the  next  chanter  must  be  reviewed. 

Furthermore,  if  the  upper  limits  of  the  buoyancy  subrange  for 

-  ^ 

the  velocity  and  temperature  spectra  are^ interest^,  the  wave  numbers 
k*  and  k*  obtained  by  equating  Eqs.  2.2  and  3.140  as  well  as  Eqs.  2 
and  3.141  respectively  are  expressed  as-  ^ ^ 


/  //m 

.3 


3(4+C4) 


3(4+C4) 

~8(2-C  )  8(2-C .)  i  -j 

k*  ~  a  f  N  e 


_S_ 


and 


3(4+C4l  3(4+C4) 

8-4C  ,  8-4C " 


/ 

\  T 

1  1 

f 

c  / 

(3.142] 


-  “T 


e 


t 


Certainly,  k*  and  k*  are  linearly  yroyortic 


(3.143) 


to  k  defined  in 
o 


pr* 


Eq.  3.97.  Hence,  it  has  been  proved  aposteriori  that  k^  is  a 


characteristic  wave  number  to  distinguish  the  buoyancy  subrange  from 
the  inertial  subrange  in  case  k^  <<  k^  ;  of  course,  if  k^  2  kQ  or 
k  >>  k  ,  there  exists  no  buoyancy  subrange.  Equations  3.142  and 


3.143  also  reveal  another  interesting  thing,  i.e.,  the  upper  limit  of 
buoyancy  subrange  are  a  function  of  flow  conditions  such  as  dT/dz^ 
etc. 
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Chapter  IV 
RESULTS  AND  ANALYSIS 


In  this  chapter,  numerical  solutions  of  turbulence  soectra  of 


a  thermally  stratified  flow  investigated  in  the  last  chapter  will  be 
given.  The  solutions  studied  in  the  previous  chapter  are  mainly 
divided  into  two  categories:  one  is  for  the  solution  of  a  flow  with 
molecular  effects,  the  other  is  for  those  of  a  flow  with  neglijgiMe 
molecular  effects.  The  consideration  of  molecular  effects  in  the 
flow  of  the  first  category  does  not  mean  that  the  spectra  at  high 
wave  numbers  beyond  the  molecular  cut-off  wave  number  can  be  studied 
by  the  present  method--the  generalized  eddy-viscosity  approximation, 
but  is  used  to  generate  nondimensionalized  spectra  expressed  with 
dimensionless  variables  containing  the  molecular  parameters  v  and 


Thus^  for  studying  huovancv  subrange  ,  we  need 

only  consider  the  case  when  the  molecular  effects  are  negligib le  in 
the  wave  number  range  of  interest.  In  other  words,  the  solution 
given  in  section  3.3.1  will  not  be  investigated  exf_an.&  i  vp  ly  and  only 
some  typical  spectral  curves  of  a  certain  flow  condition  are  displayed 
as  Figs.  28,  29,  and  30.  However,  for  a  better  understanding  of  the 
spectral  forms  in  the  buoyancy  subrange,  the  solutions  developed  in 
section*3. 3.2  must  be  ascribed  to;  in  particular  the  solution 
obtained  in  case  s  ->  °°  and  r  ->  <»  as  presented  in  part  B  of 
section  3.3.2  will  be  studied_extensively  because  of  the  clearer  and 
simpler  forms  given  b}^£gs .  3.H7  and 
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4.1  Determination  of  Parameters  r  ,  ,  r^,  ,  C.,  anJ 

In  order  to  solve  the  nondimensionalized  spectra  such  as  $(x) 

and  $^(x)  ,  etc.  as  functions  of  the  dimensionless  wave  number  x  by- 

using  Eqs .  3.117,  3.H_8,  3.119  and  3.120,  the  parameters  r  ,  T1  ,  T  , 

'  -  w  v-i /  uL 

and  must  be  known.  From  Eq.  3.96,  r  ,  r.  ,  and  T  can  be 

\  C~Yif  dU/ 


related  to  only  two  parameters  and  if  dlJ/dz  ,  dT/dz  , 

b  ,  N  ,  e  and  g  are  known/  Experimentally ,  and  can  be 

evaluated  from  the  measured  spectra  at  given  flow  conditions  charac¬ 
terized  by  dU/dz  ,  dT/dz  ,  b  ,  N  ,  e  and  g  .  Thus ,  the  way  to 
evaluate  and  would  be  equivalent  to  that  of  evaluating 

- . — \ r — -“T _ r  * 

Kolmogorov's  constant  a  in  Eq.  2.2  for  locally  isotropic  turbulent 
flow.  For  simplicity,  the  asymptotic  forms  of  Eqs.  3.132  and  3.133 


can  be  used  for  evaluating  ^  if  there  exists  a  axicy  sub¬ 
range  and  if  the  flow  conditions  upon  which  depends  are  known. 

However,  for  the  present  study,  these  parameters  can  only  be 
assumed  before  the  numerical  solutions  are  obtained.  Generally 
speaking,  dU/dz  ,  dT/dz  ,  and  g  can  be  found  from  the  measure¬ 
ments  of  velocity  and  temperature  profiles  of  the  atmosphere,  and  b  , 
N  ,  and  e  can  be  estimated  from  the  measurements  of  heat  flux  and 
shear  stress  or  from  the  energy  balance  budget  equation  or  from  the 
measured  spectra.  Thus,  we  can  evaluate  the  maximal  and  minimal 
values  of  those  flow  characteristics,  but  how  they  are  related  to 
one  another  for  given  mean  velocity  and  temperature  gradients  is 


unknown.  Nevertheless,  to  the  best  knowledge  of  the  author,  there 

Care  no  measurements  of_  N_  and e  corresponding  to  the  respective 

measurements  of  spectra  in  the  free  atmosphere.  Therefore,  in  the 
present  study,  the  values  of  r  ,  , 


and  rT  are  assumed  to 
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facilitate  the  numerical  solutions,  and  the  significance  of  these  three 
parameters  can  only  be  studied  by  varying  their  relative  values.  On 
the  other  hand,  and  are  varied  according  to  0  <  and 

0  <  in  order  to  test  how  the  spectral  forms  will  be  changed  under 


a  given  set  of  values  for  r  ,  ,  and 


Henceforth,  we  can 


obtain  a  generaj^  idea  that  the  introduction  of  and  gives 


more  degrees  of  freedom  to  test  the  spectral  forms,  although  more 

freedom  means  more  complications  are  involved  in  the  analysis.  In 

f 

the  following,  results  and  analysis  are  presented  and  classified 
according  to  the  flow  stratification. 


4.2  of  Stably  Stratified  Flow 

4.2.1  Asymg^fit^c  solutions  of  stably  stratified  flow  in  the 
buoyancy  Subrange  by  V^v-ina ,  -  As  studied  in  section  3.4,  in  the 

buoyancy  subrange  where  the  productions  of  energy  and  temperature 
inhomogeneity  are  negligible,  that  is  equivalent  to  say  that  |r|  ~  0 
and  r  >>  | r  |  in  Eqs.  3.117  and  3.118  from  the  numerical  point  of 


view,  there  exist  some  asymptotic  solutions  of  velocity  and  tempera-  ^ 
ture  spectra  as  functions  of  .  Thus,  it  is  proposed  to  solve. 

Eqs.  3.117  and  3.118  numerically  bv  letting  \v\  =0  ,  |  |  -  3^  ,  and 

|  r ^  |  be  finite^.  ("Fig^e^4^1i splays  how  the  spectral  forms  vary  as  the 
parameter  changes  for  stable  stratification.  In  the  derivation 

of  Eqs.  3.132  and  3.133,  is  restricted  to  be  0  <  ,  however, 

if  Eqs.  3.117  and  3.118  are  concerned,  needs  only  satisfy  the 

condition  <  1  when  jdT/dzj  is  very  small  as  discussed  in 

section  3.2,  and  no  lower  limit  should  be  assigned  to  from  the 

numerical  point  of  view.  Thus,  for  the  sake  of  interest,  in  Fig.  4, 
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some  curves  evaluated  by  Eqs.  3.117  and  3.118  in  case  C.  <  0  are  also 
presented.  Of  course,  the  power  law  derived  in  Eqs.  3.132  and  3.133 
cannot  be  true  in  this  case. 

As  displayed  in  Fig.  4,  the  spectral  slopes  in  the  buoyancy 

subrange  show  agreement  with  these  derived  in  Eqs.  3.132  and  3.133, 

*  „  n  nn.  .  -  11.999/4.001  -3 

i.e.,  for  =  0.001  ,  <2>(x)  ~  x  ~  x  and 

,  (  .  -  4.003/4.001  -1  ,  ^  n  n  .  -  11.1/4.9 

^^(x)  x  "  x  ,  and  =  0.9  ,  $(x)  ~  x 

and  $^(x)  ~  x  6. 7/4.9  .  jt  is  clear  that  with  =  0.001  the 

asymptotic  solution  becomes  very  close  to  the  prediction  of  Lumley- 

Shur,  and  with  =  0.9  the  slope  approaches  the  predicted  slope 

given  by  Bolgiano.  For  better  understanding  the  behavior  of  the 

spectra,  Fig^5  is  plotted  by  varying  at  given  =  0 . 00  k 

Thus,  we  can  see  that  as  T-  increases  the  buoyancy  g ffr>r.t_s.-p^neJ:rate 

“  uLr 

gradually  into  higher  wave  numbers  and  the  deviations  from  the  -  5/3 

■■ »  /—  ■  ■  '  i  mm m- 

law  of  locally  isotropic  flow  are  more  apparent  as  the  spectra  become 

more  anisotropic  due  to  the  buoyancy  effect s  character i r ^  . 
Figure  31  displays  the  asymptotic  power  law  varied  as  for  velocity 

and  temperature  spectra  in  the  buoyancy  subrange  of  stably  stratified 
flow,  m  and  n  contained  in  c{> (k)  -  kn  and  ^^(k)  -  km  are 


K* 


defined  as  m  = 


-3C.  -4 
4 

4+C 


n  = 


C  -12 
4 

4+C , 


4  4 

4.2.2  The  production  effects  of  turbulence  energy  and 
temperature  ^ixi^ppgcqeity^  on  the  spectra  of  stably  stratified  flow 
in  the  buoyancy  subrange  -  In  section  4.2.1,  it  was  assumed  that  the 
production  of  turbulence  energy  and  temperature  infromogeneity  is 
negligible  in  the  buoyancy  subrange.  Without  any  doubt,  in  this 
buoyancy  subrange,  only  the  inertial  transfer  process  and  the  buoyancy 
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f 


Ui 


k 


effects  exist  and  the  spectral  forms  must  characterize  the  buoyancy 
effects.  However,  if  |r|  and  |r  |  are  large  enough  as  compared 
to  I\  so  that  the  production  effects  can  penetrate  into  the  buoyancy 

n _ is - -  - =■  " 

subrange ,  the  spectral  forms  will  hp  Hi stnrheH .  Figure  _6  shows  two 
sets  of  spectral  curves  at  different  r^,  ,  r  ,  and  with  fixed 

C,  and  .  The  (J^^J^e^jDf^curve^  i-s  calculated  with^^^=  0 . 0 1  , 
^£^=^0^  ,  fj  =  0.1  ,  L'2  =  1.0  ,  and  and  the  with 

(kT  =  0.001  ,  ^  =  0.001.  r,  =  0.01  t>  C?  =  1.0J  ,  and  C4  =  0 .jj  From 
the  plotted  curves,  the  temperature  spectra  present  maximal  points 
which  show  that  some  energy  is  converged  fronL the  yplncity.  fiejd  into 

the  tem^£rature^frekl .  Nogr  wav^  nimihp.r  r;lPfp  in  which  the  maxi¬ 
mum  of  tomperatur^—snectr^  occurs,  there  appears  a  steeper  slope  than 
-^^^^or^the  energy  spectrum  as  a  consequence  of  kinetic  energy^  being 
transferred  and  converted ^intoS pptentjal  energy^  It  is  clear  that 
this  wave  number  range  is  associated  with  the  so-called  buoyancy  sub¬ 

range  in  our  previous  investigation  of  asymptotic  solutions. 

According  to  the  asymptotic  solutions  presented  in  section 
3. 4.  L»  the  energy  spectrum  and  thq  temperature  spectrum  must  show  the 

~&¥  C^5>  — 

power  law  -  11.7/4.3,  and  -4. 9/4. 3  ^respectively  if  only  the  buoyancy 
effects  are  predominant  in  the  buoyancy  subrange^  however,  the  pqwer 
law  of  velocity  spectrum  of  the  first  set  of  curves  appearing  in 
. 4  .  This  is  not  surprising  after/^igT  ^^is  reviewed. 

Figure  7  displays  the  distributions  of  energy  production,  transfer  and 
drainage  by  buoyancy  force  as  represented  by  ,  F0  ,  and  F^ 


jP£o2>. 


(4.1) 


T£fiUSF. 


(4.2) 


5  3 


pMWIlCF , 


(4.3) 


The  functions  and  are  the  production  and  inertial  transfer 

of  temperature  inhomogeneity,  figure  7^hows  that  in  the  wave  number 
range  the  -  5/3  law  holds  the  inertial  transfer  presented  as  F^ 
predominates  the  turbulent  structure  and  tends  to  be  a  constant.  In 
case  the  energy  drainage  by  vertical  heat  flux  becomes  gradually  im¬ 
portant,  F0  increases  as  wave  number  decrease^  and  consequently  the 
buoyancy  subrange  is  formed.  As  the  wave  number  decreases  down  to 
the  region  where  the  energy  production  becomes  important 
decreases.  Now,  in  the  wave  number  range,  say  in  the  interval 

* 

0.03  <  x  0.15  ,  the  energy  draina£Q_by,  buoyancy  becomes  important 

and  the  energy  production  becomes  less  negligible  as  well.  In  other 
words,  the  introduction  of  the  energy  production  in  this  wave  number 
range  compensates  the  energy  drainage  by  buoyancy  and  modifies  the 

power  1 aw  f r om  ^ -  .11. 2JA  * 3  La _ - 2 . A  .  CA  similar  situation  can  be 

seen  from  the  second  set  of  spectral  curves  which  shows^  that  the 
reduction  in  energy  production  can  cause  a  wider  buoyancy  subrange 


and  steeper  power  law.  for  a  better  understanding.  Fig.  8  must  be 

— - - ^ - - - - J 

examined.  In  the  interval  0.004  <_  x  <_  0.02  ,  the  turbulent  energy 

is  mainly  distributed  by  energy  transfer  F^  and  drained  by  buoyancy 
F^  but  the  contribution  from  energy  production  F^  on  the  energy 
distribution  is  less  important.  Hence,  it  will  not  be  a  surprise  to 
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have  the  maximum  power  law  -2.5  for  the  energy  spectrum  curve  B  in 
the  wave  number  range  0.004  <  x  <  0.02  . 


lother  interesting  phenomenon  also  caused  by  the  compensation 
of  energy  product i on ^on  F0  and  can  be  observed  in  the  wave 
number  range  x  <  0.001  where  the  -1  law  exists  for  the  curve  _A  and 


the  -5/3  law  for  the  curves  B  in  Fig.  6.  /No>^referring  tc^Fig.  7^ 
F^  completely  predominates  the  energy  distribution  in  wave  number 
range  x  <  0,001  ,  the  -1  slope  for  jurve  A  \can  be  predicted  from 


Ec|.  4.1  by  inserting  *  As  to  v^rve  Bjt  -5/3  slope  in  the  wave 

number  x  <  0.001  results  from  the  condensation  among  F,  ,  F^  , 

■  1  i  1.  i  4 


and  F^  such  that  F^  tends  to  be  constant  over  a  wide  range  as  can 


be  seen  from  Fig^_  8.  Here  one  word  must  be  added,  the  -5/3  slope  in 
the  wave  number  range  x  <  0.001  is  obviously  not  a  result  of  local 


isotropy.  Thus,  to  the  experimenters,  the  prediction  for  the  local 
isotropy  from  the  measured  velocity  spectra  must  be  carefully  worked 
out  in  case  buoyancy  effects  exist.  For  example,  if  in  the  buoyancy 
subrange  there  exists  any  experimental  error  which  causes  data 
scattering,  the  -5/3  slope  may  be  extended  to  low  wave  number  without 
realizing  the  existence  of  the  buoyancy  subrange.  In  other  words, 
the  -5/3  slope  appeared  in  low  wave  number  x  <  0.001  for  curve  B 
in  Fig.  6  may  incorrectedly  be  predicted  as  a  result  of  local  isotropy. 
The  experimental  data  of  Fig.  2a  shows  this  situation  very  clearly 
since  a  -5/3  slope  appears  on  the  left  side  of  the  buoyancy  subrange. 

In  order  to  test  the  local  isotropy  of  the  stably  stratified 

5/3 

turbulent  shear  flow,  Fifl ._9  is  plotted  with  x  $(xj)  and 
x5^3^tt(x)^  vs  x  It  is  clear  that  the  plot  of  Fig.  9  can  present 

the  tendency  to  local  isotropy  in  a  better  way  than  the  plots  0(x) 
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and  $TT(x)  vs  x  ,  since  the  deviation  from  local  isotropy  can  be 
detected  easily  without  any  ambiguity  as  stated  in  the  last  paragraph. 

4.2.3  Validity  of  the  generalized  eddy-viscosity  approximation- 
Up  to  this  stage,  the  generalized  eddy-viscosity  approximation  has  been 
used  to  study  turbulent  energy  and  temperature  spectra  under  the 
effects  of  buoyancy.  Of  course,  the  validity  of  the  present  study 
can  be  justified  from  the  measured  energy  and  temperature  spectra, 
however,  there  is  another  way  to  test  its  validity  by  comparing  the 
shear  stress  and  vertical  heat  flux  spectra  predicted  by  the  present 
study  with  respect  to  the  measured  ones.  Since  the  present  method-- 
generalized  eddy-viscosity  approximation--rests  on  the  assumptions 
for  the  integrated  forms  of  energy  transfer,  temperature  inhomogeneity 
transfer,  shear  stress,  and  vertical  heat  flux  spectra  presented  as 
Eqs.  3.56,  3.58,  3.61,  and  3.62,  the  validity  of  the  present  study 
must  logically  be  checked  by  comparing  those  spectra  with  the  measured 
ones  although  the  measurements  of  these  spectra  are  not  easily 
performed . 

Thus,  for  the  above  reasons,  some  spectra  of  shear  stress  and 
vertical  heat  flux  predicted  by  the  present  stud)"  for  given  flow  con¬ 
ditions  are  displayed  in  figs.  10  And^l  1 ,  In  Fig.  11,  the  heat  flux 
spectrum  4^(x)  of  curve  A  shows  a  change^in  _sign  at  wave  number 
x  =  0. 00315,.  The  same  phenomenon  also  occurs  in  Fig.  12  which  shows 
the  Bolgiano's  solution  in  case  =  1  .  From  these  results  the 
changing  sign  of  $^(x)  seems  to  be  related  to  . 

For  further  understanding  turbulence  spectra  of  stably 
stratified  flow ,  Figs .  l^JiL  are  displayed  systematically  by  varying 
one  of  the  parameters  ,  T  ,  when  the  others  are 
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fixed.  In  the  spectra  the  buoyancy  subrange  certainly  exists  since 

for  most  the  condition  r  >>  r  is  satisfied.  For  example,  in 
_  w  I 


Fig.  13* *  when  is  increased  from  0.001  for  curve  A  to  0.01  for 

curve  B  with  the  corresponding  =  0.5  ,  the  buoyancy  subrange  is 
narrowed  as  expected. 

Fig.  14»  the  effects  of  energy  production  on  the  spectra 
are  examined  by  varying  r  and  keeping  the  other  parameters  fixed. 
Increasing  T  not  only  narrows  the  buoyancy  subrange  but  even  shrinks 
the  inertial  subrange.  This  situation  can  be  seen  from  the  curve  A 


of  velocity  spectra.  Near  the  region  x  =  1  ,  the  slope  of  velocity 
spectrum  is  -5/3  ,  but  in  the  interval  0.1  <  x  <  1  ,  the  slope  is 
l£ss  than  -5/3  ,  because  the  effects  of  energy  production  penetrate 
deeply  into  the  region  of  high  wave  numbers  such  that  the  energy 
transfer  decreases  with  decreasing  wave  number  in  this  interval. 
While  in  the  buoyancy  subrange  of  stably  stratified  flow  the  energy 
transfer  increases  with  decreasing  wave  number  and  the  slope  of  the 
velocity  spectrum  is  greater  than  the  slope  -5/3  .  This  situation 
reflects  the  case  when  the  flow  has  great  shear  gradients 


Figure  15>displays  a  case  when  T ^  varies.  From  the  plot, 
increasing  T  would  mean  that  the  effects  of  buoyancy  force  inten- 

*v! 

sify  as  indicated  by  a  wider  buoyancy  subrange  shown  as  curve  C. 

Also  as  r  increases,  the  power  law  in  the  buoyancy  subrange 
approaches  to  the  asymptotic  power  law  (C4-12) /(4+C4)  .  Another 


interesting  thing  should  be  noted  from  the  plots  For  these 

spectra,  C0  is  varied  and  we  can  see  that  when  ,  the 

buoyancy  subrange  for  the  velocity  spectrum  <HX)  disappears  as  C ^ 
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decreases  to  0.1  for  curve  A.  Whether  this  phenomenon  is  realistic 
or  not  can  only  be  determined  by  experiment. 

fn  Fig.  f7J  is  varied.  As  we  can  see  that  as  varies 

from  0.5  for  curve  A  to  1.0  for  curve  B,  the  Jiump  in  $^(x)  dis¬ 
appears.  Also  since  T  is  not  negligible  as  compared  to  T and 
T ^  ,  the  power  law  in  the  buoyancy  subrange  deviates  greatly  from  the 
asymptotic  power  law  which  requires  that  T  is  negligible.  In 

is  forced  to  be  negative  and  a  steeper  slope  of  the 
velocity  spectra  in  the  buoyancy  subrange  is  observed.  In  case 

=  -  0.1  the  slope  is  close  to  the  power  law  for  $(x)  given  by 
(C^-12) / (4+C^)  ,  but  for  =  -  0.5  the  plotted  curve  has  a  steeper 
slope  than  (C^-12)/(4+C^) .  This  is  not  unexpected  since  the  power 
law  derived  in  bq.  3.125  is  valid  only  for  >  0  . 

4 . 3  Buoyancy  Subrange  of  JJastaM^^t^atification 

As  described  in  Eqs .  3.127  and  3.128,  the  buoyancy  subrange 
o f  fins t ab lyj  s t rat i f ied  turbulent  flow  exists  when  the  product ion  of 
turbulent  energy  is  less  important >>  |r^|  ,  and  the  velocity 
spectrum  exhibits  a  hump  and  the  temperature  spectrum  has  a  steeper 
slope  than  -5/3  in  the  buoyancy  subrange.  Figure (T^di splays  spectral 
curves  for  unstably  stratified  flow.  Curves  A,  B  and  C  show  that 
buoyancy  effects  very  clearly  since  T ^  >>  |r^|  is  certainly  ful¬ 
filled,  however,  curves  D  appear  in  a  different  way  from  curves  A,  B 
and  C  just  because  of  the  fact  that  =  r  =  0.001  . 

In  Fig.  the  effects  of  the  production  of  turbulent  energy 

due  to  Reynolds  stress  on  the  spectra  are  displayed  by  varying  T 
when  the  other  parameters  are  fixed.  As  can  be  seen  from  Fig.  20, 
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the  variation  of  T  does  not  affect  the  spectra  significantly, 
however,  if  is  varied,  i.e.,  the  degrees  of  interaction  between 

mean  strain  rate  and  the  eddy  strain  rate  are  varied,  the  spectra 
appear  differently  as  presented  by  Fig.fj^>  A  plot  of  curves 

(i  1,...,5)  vs  x  is  shown  in  Fi^^2^7when  F^  represents  the 
terms  on  the  right  sides  of  Eqs .  3.117  and  3.118  respectively.  F^ 
decreases  with  decreasing  wave  number  as  a  characteristic  of  buoy¬ 
ancy  subrange  of  an  unstably  stratified  flow.  Note  that  in  the  bypy- 
ancy  subrange  of  a  stably  stratified  flow.  F0  the  energy__transfer 
flux  increases  with  decreasing  wave  number  as  a  consequence  of  com- 
pensating  the  energy  drainage  by  vertical  heat  For  a  better 

understanding  of  how  the  shear  stress  spectrum  is  distributed  with 
respect  to  wave  number  in  these  cases  of  unstably  stratified  flow, 
Fig.  ,2$  is  displayed  by  plotting  x^^(x)  vs  x  on  semi logri thmic 
paper.  The  interesting  thing  is  the  changing  sign  of  >  °f 

course,  the  validity  of  the  curves  of  x^^(x)  plotted  in  Fig.  23 
must  be  checked  from  experiments  although  it_  is  certain  that  ^uw(x) 


need  not  necessarily  be  always  positive  or  negative  throughout  all 
wave  numbers.  In  addition,  Fig.  (24^)is  presented  for  the  heat  flux 
spectrum  which  does  not  show  any  change  in  sign. 

The  curves  discussed  above  mainly  correspond  to  the  flow 
with  Tl  >_  |  rT|  .  In  Fig  the  case  |r  |  >>  r  is  given  by 

varying  C0  and  is  used  for  comparison  with  Fig.  21.  Up  to  now, 
the  plotted  curves  of  temperature  spectrum  seem  aJLways  to  have 
steeper  slope  than  -5/3  ,  however,  Fifl.  ffrUhows  that  this  may  not 
always  be  the  case  when  some  special  values  of  the  parameters  are 

considered;  of  course,  whether  the  curves  shown  in  Fig.  26  are 
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realistic  or  not,  or  say,  whether  the  assigned  values  of  parameters 
are  reasonable  or  not  must  be  determined  by  experiment. 

In  sections  4.2  and  4.3,  great  attention  has  been  paid  to  the 
characterist ics  of  the  buoyancy  subrange  of  both  stable  and  unstable 
stratifications.  However,  very  often  there  exists  no  buoyancy  sub¬ 
range  at  all  because  of  the  fact  that  the  criterion  F^  >>  |  F^  |  for 
buoyancy  subrange  is  not  satisfied  and  that  the  effects  of  the  pro¬ 
duction  of  turbulent  energy  penetrate  greatly  into  high  wave  number 

range^  For  example,  in  case  =  1  ,  C-  =  1  ,  the  -1  law  exists  for 

~~  - - xl - -J  L 1 _ J  ~  - - 

both  velocity  and  temperature  spectra  as  the  productions  of  energy 

and  temperature  inhomogeneity  are  very  large.  Figured^  shows  this 

situation  for  both  stratifications.  Also  in  Fig.  27,  the  curves  A 

show  the  existence  of  - _l/3  slope  region  for  temperature  spectrum 

when  the  velocity  spectrum  has  extensive^-  5/3  slope  in  this  region. 

This  region  may  be  called  the  inertial  and  convective  region  as 

stated  in  Gisinafs  paper  ( 1966)  and  as  indicated  by  Eq.  3.129  when 

C  .1  ' 

^4  • 
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Now  in  order  to  examine  systematically  how  the  parameters  F  , 
F^  ,  F^  ,  C 2  ,  ,  s  ,  and  r  vary  for  these  figures,  a  brief  table 
is  listed  as  follows: 
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Now,  the  general  characteristics  of  turbulence  spectra  of 
thermally  stratified  flows  has  been  investigated  very  extensively  for 
the  buoyancy  subrange  by  means  of  the  generalized  eddy-viscosity 
approximation,  it  would  be  worthwhile  to  make  comparisons  with  mea¬ 
surements  in  the  atmosphere  and  to  test  the  validity  of  the  general¬ 
ized  eddy-viscosity  approximation. 


4.4  Qualitative  Cgpggja^n  Between  the  M,e^snr^matt4s  of  Spectra  in 

the  Atmosphere  and  the  Results  Obtained  by  the  Generalized  Eddy- 

Basically  speaking,  it  is  hard  to  make  this  comparison  because 
of  the  fact  that  the  experimental  data  normally  relate  to  one- 
dimensional  spectra  whereas  the  theoretical  consideration  in  the 
present  study  is  concerned  with  three-dimensional  spectra.  As  noted 
in  the  paper  of  Alkseev  and  Yaglom  (1967) ,  the  one-dimensional-jsp^c- 
trum  always  varies  noticeably  morg  smoothly  than  the  three- 
dimensional  spectrum.  The  wave  number  range  in  which  a  power  law 
occurs  in  a  one-dimensional  spectrum  will  not  be  the  same  range  in 
which  the  same  power  law  occurs  in  a  three-dimensional  spectrum  even 
if  the  local  isotropy  of  the  turbulence  holds.  l>1o^ove^^  the  turbu- 
lence  spectra  of  a  thermally  stratified  £low^are  certainly  charac¬ 
terized  b  s  to  the  buoyancy  effects ,  and  thus  in 

our  nroblem.  the  basic  advantage  of  local  isotropy  to  relate  one- 
dimensional  spectra _to three-dimensional  spectra  is 

However,  if  only  the  approximat  ^comparison  is  made  to  see 
how  the  parameters  mentioned  in  the  theoretical  consideration  are 
related  to  the  experimental  data,  the  difficulties  as  stated  above 


can  be  relaxed.  In  particular,  it  is  assumed  that  in  the  buoyancy 
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subrange  the  deviation  from  the  local  isotropy  may  not  be  too  great. 

In  section  2.4,  it  is  mentioned  that  the  exponent  n  of  the  velocity 
spectrum  k  n  in  the  buoyancy  subrange  of  a  stably  stratified  flow 
varies  from  2.0  to  3.5.  cAs  we  reviev^Fig.  4y  the  values  of  n  are 
associated  with  the  parameter  £k)  in  the  theoretical  consideration, 
i.e.,  by  varying  ,  some  specific  s lope  in  the  buoyancy  subrange 

can  be  made  or  example,  to  n  _=  liys  there  corresponds  =  1  , 
for  n  =  3  we  have  -*  0  ,  and  for  higher  values  of  n  ,  a  negative 
value  can  be  assigned  to  as  presented  in  Figs.  4  and  18.  And  in 

fact,  the  introduction  of  C.  into  the  theoretical  consideration  is 
a  great  improvement  in  interpreting  the  spectra  of  a  stably  stratified 
flow. 


At  tjxis..  st.ag&^  it  would  be  worth  mentioning  jthe  process  by 
which  the  parameter  was  introduced  into  the  theoretical  considera¬ 


tion.  In  the  inertial  subrange  of  velocity  spectrum,  only  e  the 
dissipation  of  turbulent  energy  is  the  parameter  to  characterize  tur¬ 
bulence.  Thus,  the  -  5/3  law  is  implied  from  the  dimensional  argument 
and  this  subrange  is  of  universal  equilibrium  since  the  flow  is  inde¬ 
pendent  of  the  parameters  such  as  mean  velocity  and  temperature  gra¬ 
dients  which  in  turn  characterize  the  external  flow  conditions.  As 
to  the  energy  containing  range,  the  flow  is  determined  by  the  mean 
quantities  and  is  influenced  by  the  geometry  which  contains  the  flow 
or  around  which  the  flow  passes.  Of  course,  this  range  cannot  be 
universally  determined,  and  moreover,  there  exists  no  equilibrium 
state  for  turbulent  flow  in  this  range.  Now,  if  we  assume  that  there 
exists  a  certain  subrange  between  the  inertial  subrange  and  the  energy- 
containing  range,  then  we  can  expect  that  the  turbulence  spectra  are 
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in  a  state  of  equilibrium  but  toot)  of  universal  form.  Consequently, 


it  can  be  expected  that  in  this  subrange  the  turbulence  is  determined 
by  the  mean  quantities  and  the  turbulence  parameters. 

^First,  consider  a  subrange  which  is  caused  by  the  velocity 
field  alone  such  as  mean  velocity.  Certainly  dU/dz  and  e  are 
the  only  parameters  to  determine  the  turbulence  as  can  be  seen  from 
Tchen's  solution  (1953).  ^Seco^id^  if  this  subrange  is  induced  by  the 
existence  of  thermal  stratification,  then  based  on  the  previous  argu¬ 
ments,  the  turbulence  in  this  "buoyancy  subrange”  must  be  determined 
by  the  mean  quantities  dT/dz  ,  g/T  ,  and  the  turbulent  parameters 
e  and  N  .  Note  in  the  present  arguments  the  significance  of  the 
second  and  higher  order  derivatives  of  U  and  T  with  respect  to  z 
are  negligible  as  compared  to  dU/dz  and  dT/dz  respectively.  The 
following  paragraph  gives  the  reason  why  the  parameter  C^  must  be 
introduced. 


In  the  buoyancy  subrange  the  turbulence  may  be  in  equilibrium 
but  not  universal  in  character.  Thus ,  for  a  given  set  of  parameters, 
dT/dz  ,  g/T  ,  e  ,  and  N  of  a  certain  turbulent  flow,  there  exists 
a  definite  spectral  form  beqause  of  the  equilibrium  of  the  turbulence. 
However,  the  spectral  form  will  vary  from  one  turbulent  flow  to  the 
other  because  the  turbulence  is  not  universal.  Hence,  we  can  see  the 
necessity  to  introduce  a  new  dimensionless  parameter  /~C^j  in  order 
to  characterize  the  spectral  forms  for  varied  flow  conditions  as 
related  to  the  degree  of  interaction  between  the  mean  temperature 


field  and  temperature  fluctuation  field  as  described  in  Eq.  3.62 .  Of 


course,  the  same  argument  can  be  applied  to  )  C 


introduced  in  Eq.  3.61. 


In  >  TchenVs  solutions  (1953)  are  only  two  particular  cases  implied 
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w  isjcnown  for  a  given  set  of  flow  conditions,  is  it 
possible-t^-f-ind  a  similarity  theory  such  that  some  characteristic 
variables  can  form  a  set  of  dimensionless  variables  and  the  turbu¬ 


lence  becomes  qu  a,s  i^un-  i  vers  a  1  in  this  subrange?  Here  the  term  "quasi- 
universal"  means  that  for  different  flow  conditions  which  can  build 
up  the  same  values  of  dimensionless  parameters  there  is  a  unique 
spectral  form  corresponding  to  the  associated  .  Certainly,  the 

characteristic^_Lo  be  quasi-universal  in  the  buoyancy  subrange  of  a 
stahl^-s^atified  flow  was  ignored  and  undetected  by  Bolgiano  and 
Lumley  since  both  of  them  devoted  themselves  to  an  effort  to  reach  a 
universal  solution  for  the_  buoyancy  subrange.  Thus,  the  reason  why 
some  dimensionless  parameters  are  presented  in  Eqs.  3.65,  3.74,  3.96, 
3.97  and  3.100  can  be  seen  clearly.  With  these  ideas  in  mind,  it 
should  i KTL  he  -surprising  to  get  the  hypothesis  (2)  for  the  modified 
hypotheses  stajLed— Ln_ s ecfipn  3.5. 

Now,  it  is  time  to  clarify  why  the  parameter  dlJ/dz  does  not 
appear  in  the  hypothesis  (2).  In  Monin’s  paper  (1965)_  it  was  noted 
that  in  the  buoyancy  subrange  the  dependence  of  spectra  upon  dT/dz 


only  and  hot  upon  dU/dz  seems^nnatural.  However,  if  dlT/dz  is 
introduced  into  hypothesis  (2)  in  order  to  make  this  hypothesis  com¬ 
plete,  the  simple  solution  as  expressed  in  Eqs.  3.132  and  3.133 
cannot  be  obtained.  However,  the  situation  will  not  be  so  pessimis¬ 
tic  since  the  dependence  of  the  spectra  upon  dU/dz  can  be  investi¬ 
gated  numerically,  although  not  analytically.  And  in  fact,  the 
numerical  investigation  of  the  dependence  of  spectra  upon  dU/dz  is 
certainly  an  improvement  upon  either  Bolgiano’s  or  Lumley-Shur f s 

- ■  |  Will  ■  I  n  |  ,,  — !!■  ™  ■  ■ 

hypothesis.  Curve  A  in  Fig.  6  shows  the  effect  of  dU/dz  on  the 
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spectra  clearly.  If  Fig.  2b  is  reviewed,  the  less  steep  part  of  the 
curve  on  the  left  side  of  the  buoyancy  subrange  reflects  the  effects 
of  the  production  of  turbulent  energy  due  to  Reynolds  stress  or  flux 
divergence. 

In  case  the  unstable  stratification  is  concerned,  the  velocity 

<  m 

spectra  shown  in  Fig.  3  show  a  hump  at  the  frequency  corresponding 
to  the  motion  of  Benard  cell  in  the  thermal  convection  (Ivanov  and 
Ordanovich,  1967),  In  the  present  study,  this  kind  of  hump  due  to 
unstable  stratification  can  also  be  detected  from  Fig.  19,  and  it 
seems  to  the  author  that  the  left  side  of  the  hump  of  the  spectra  in 
Fig.  3  approaches  to  the  +1  slope  as  shown  in  Fig.  19  from  the  theoret¬ 
ical  consideration.  Unfortunately,  Ivanov  and  Ordanovich  did  not 
mention  any  measurements  of  the  temperature  spectra  corresponding  to 
the  flow  conditions  under  which  the  velocity  spectra  were  taken. 

Hence,  generally  speaking,  when  compared  to  the  above  stated  measure¬ 
ments  of  spectra  in  the  atmosphere,  the  theoretical  study  by  the 
present  method--generalized  eddy-viscosity  approximation--can  give 
better  features  of  turbulence  structure  of  thermally  stratified  flows 
than  any  of  the  previous  hypotheses  of  Bolgiano  and  Lumley  in  the 
case  of  stable  stratification  and  than  Monin fs  results  (1962)  on  the 
spectra  of  unstably  stratified  flow. 
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Chapter  V 

CONCLUSIONSjXND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

In  this  theoretical  investigation,  it  was..sWn  f hat  the 
proposed  generalized  eddy-viscosity  approximation  can  predict  features 
of  turbulence  spectra  of  thermally  stratified  turbulent  shear  flows 
better  than  any  previous  hypotheses,.  The  present  method  gives  a 
general  solution  in  the  buoyancy  subrange,  and  the  effects  of  the 
turbulent  energy  production  in  the  buoyancy  subrange  are  investigated 
numerically.  As  a  consequence  of  the  application  of  the  generalized 
eddy-viscosity  approximation,  modified  hypotheses  arg  ahl  i  <;hed  fpr 
the  buoyancy  subrange  of  a  stably  stratified  flow  as  follows:^ 

The  buoyancy  subrange  of  a  stably  stratified  turbulent 
flow  exists  when  the  local  production  and  local  dissipation  of  turbu¬ 
lent  energy  are  not  important  and  when  the  flow  conditions  satisfy 
the  criterion 


I  E  ldTl 

(  '  dz 

-  <<  1 

N 


g  i  ciT  j  N 

or  —2—  !  — |  <<  - 

~  dz  e 


(-H 


1/2 


which  means  that  the  internal  frequency  (N/e)  g/T .  is  much  greater 

,  dj  .1/2 

than  the  Brunt-Vaisala  frequency,  1-5- 

(2)  In  the  buoyancy  subrange,  the  spectra  of  turbulent  energy 


and  temperature  fluctuation  are  determined  hv  N  ,  e  ,  g/T  ,  dT/dz  , 


k  ,  and  completely,  and  are  expressed  in  forms 

F  -  <*— -  8_4C 
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-1+C, 


1-C, 


r,  =  n_ 


T  / 


-1+C,  1-C, 

4  AT  4 

I— I 

'dz^  1 


C,  >  0 
4 


o°  >  s  >  2 


N  =  2N  £_  =  ?  £L 

-*  ’  Jz_  dz 


and  f  and  f?  are  numerical  variab les  .  Clearly,  both  the  solution 
of  Bolgiano  and  the  solution  of  Lumley-Shur  are  contained  in  these 
forms  when  =  1  and  0  respectively. 


As  to  the  unstably  stratified  flow,  the  velocity  spectrum 
exhibits  hump  in  the  buoyancy  subrange.  On  the  left  side  of  this 
hump  the  velocity  spectrum  approaches  a  +1  slope  and  the  temperature 
spectrum  shows  a  -_5  slope. 

- —  ~  Before  the  generalized  eddy-viscosity  approximation  is 
introduced,  the  basic  assumption  used  to  derive  a  simpler  set  of 
spectral  equations  is  that  the  flow  is  locally  homogeneous.  In  other 
words,  in  the  derivation  of  the  spectral  equations  the  terms  due  to. 
inhomogenei tv  of  the_flow  field Jij^e  been  discarded.  Thus,  for  better 
understanding  of  the  effects  due  to  the  inhomogeneity^ it  would  be 
constructive  to  include  these  terms  in  the  spectral  equations.  This 
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may  be  important  in  the  case  of  free  convection  and  inversion  layer 
flows  because  the  divergences  of  energy  flux  and  temperature  inhomo¬ 
geneity  play  an  important  role  in  the  balance  of  the  thermal  turbu¬ 
lence  energy  budget  equation.  Thus,  for  further  research  on  turbulence 
spectra  of  thermally  stratified  flows,  these  effects  must  be  considcxad-^ 
so  that  the  spectral  forms  in  the  subrange  where  the  production  and 
diffusion  of  turbulent  energy  play  the  principal  role  in  determining 
the  turbulent  structure  can  be  investigated. 

On  the  other  hand,  the  fruitful  results  obtained  up  to  now 
may  induce some  further  advanced  research  in  applied  engineering 
problems  such  as  turbulent  diffusion  related  to  air  pollution, 
wave  propagation  in  the  atmosphere,  and  high  speed  aircraft  design. 

For  example,  in  his  paper,  Tchen  ( 1959)  did  jnot  consider  -the  effects 
of  stratifi cation  on  the  di s persion  of  smok^  froiq  a  j3oi n t  source 
although  the  role  of  the  Reynolds  stresses  in  the  momentum  equations 
was  studied  extensively  in  order  to  clarify  the  internal  interaction 
of  the  diffusing  particles.  In  the  classical  theory  of_  the  scattering. 
of  sound  waves  (Tatarski,  1961),  the  vertical  heat  flux  spectrum  is 
assumed  to  be  zero  in  order  to  simplify  the  problem;  certainly  the 
vertical  heat  flux  spectrum  is  not  zero  in  the  case  of  stratified 
flows.  As  to  the  high  speed  aircraft  design,  the  airplane  frequency- 
response  function  can  predict  the  associated  response  spectrum  of 
airplane  vibration  caused  by  turbulence  (Steiner  and  Pratt ,  1967) . 
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(Pinus  and  Shcherbakova,  1966) 
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Fii>.  2  Typical  spectral  density  curves  of  the  horizontal  wind  velocity 
component  of  stably  stratified  jQows  (a)  Vinnichenko,  1966 
(b)  Pinus ,  196^5^ — - 
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Fim.  3  Typical  spectral  densities  of  horizontal  wind  velocity 

components  of  unstably  stratified  flow  at  various  heights 
(Ivanov  and  Ordanovich,  1967) 
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Fig.  4  Asymptotic  velocity  and  temperature  spectra  of  stably 
stratified  flows  varied  as  a  function 
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Fig.  5  Asymptotic  velocity  and  tempexature  spectra  of  stably 
stratified  flows  varied  as  a  function  of  T ^ 
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Fig. 


6  and  temperature  spectra  of  stably  stratified 

flows  with  the  effects  of  the  production  of  turbulence 
'energy  and  temperature  inhomoj»eneit>F“  ' 
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Fig.  10  Turbulence  spectra  of  shear  stress  and  vertical 
heat  flux  of  stably  stratified  flows 


(x)  ,-<X>wT(x) 
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Fig.  11  Turbulence  spectra  of  shear  stress  and  vertical 
heat  flux  of  stably  stratified  flows 


<X>(x),  <X>TT(x  )  ,  -<tow(x)  ,  -<I>wT  C  x  ) 
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Turbulence  spectra  4>(x),  (x),  4>  (x) , 

4>wT(x)  vs  x  1  uw 


Fig.  12 


<X>(x)  i  <I>TT(x) 


in 


Fig.  13  Velocity  and  temperature  spectra  of  stably 
stratified  flows  varied  as  a  function  of 


<X>(x)  # 
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Fig.  14  Velocity  and  temperature  spectra  of  sd-ah-Ly 
stratified  flows  varied  as  a  function  of  r 


(*)  ii<j>  ‘  (X)4> 
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Fig.  15  Velocity  and  temperature  spectra  of  stably 
stratified  flows  varied  as  a  function  of 


<J>  (x  )  ,  <I>TT  (x  ) 
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Fig.  16  Velocity  and  temperature  spectra  of  stably 

stratified  flows  varied  as  a  functioiTof  C9 


<X>  (x )  , 
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Fig.  17  Velocity  and  temperature  spectra  of  stably 

stratified  flows  varied  as  a  function  of  C. 
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Fig.  18 


Velocity  and  temperature  spectra  of  stably 

stratified  flows  varied  as  a  function  of  C. 
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(x  ),  <X>  (x) 
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Fig.  19  Velocity  and  temperature  spectra  of  ups t ably 
stratified  flows  varied  as  a  function*" of  r 


U>  ,  <J>  U) 
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Fig.  20  Velocity  and  temperature  spectra  of  unstably 
stratified  flows  varied  as  a  function  of  F 


(x),  *<x) 
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Fig.  21 


Velocity  and  temperature  spectra  of  unstably 
stratified  flows  varied  as  a  function  of 
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rig.  22  Distributions  of  energy  production,  transfer,  and  input  by  buoyancy  force,  an. 

temperature  inhomogeneity  production  and  transfer  of  unstably  stratified  flow 
across  wave  numbers 


<$uw(x  ) 
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Fig.  23  x  4>uw(x)  vs  x  with  varied  C,  for  unstably 
stratified  flows 
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Fig.  24  Turbulence  spectrum  of  vertical  heat  flux  of 
unstably  stratified  flow  varied  as  a  function 
of  C2 


(x)  ,  <D(x) 
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Fig.  25  Velocity  and  temperature  spectra  of  unstably 
stratified  flows  varied  as  a  function  or""t„ 
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Fig.  26  Velocity  and  temperature  spectra  of  unstably  stratified  flows  varied  as 
function  of  CL 


4>U)  , 
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Velocity  and  temperature  spectra  of 
stratified  flows  varied  as  a  function"  of 


tably 


Fig.  27 
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Fig.  28  Turbulence  spectra  of  a  stably  stratified 
viscous  shear  flow  "" 
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Fig.  29  Distributions  of  energy  dissipation  ,  transfer  F.,  ,  drainage  by  buoyancy 
and  production  F^  ,  and  temperature  inhomogeneity  dissipation  F^  ,  transfer 
production  F  ,  of  a  stably  stratified  viscous  shear  flow  across  wave  number 


4>(x),  <*>TT(x ) 
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Fig.  30  Velocity  and  temperature  spectra  of  thermally 
stratified  viscous  shear  flows 
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